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SUMMARY 


The requirement of a low frontal area for a ram-jet has made 
necessary air velocities sometimes greater than 300 ft. per sec. 
entering the ram-jet combustion chamber. These high velocities 
make combustion difficult, and a considerable amount of research 
has been, and still is, necessary to establish the design principles 
involved in building satisfactory ram-jet combustors. The prob 
lems of combustor development can be divided into choice of the 
best fuel, proper preparation of the fuel-air mixture, ignition of the 
mixture, and completion of combustion. Factors of interest in 
choosing the fuel are ease and safety of handling, availability, 
cost, ease of combustion, volatility, and heat of combustion per 
pound and per unit volume. Data are presented showing that 
petroleum hydrocarbons are more difficult to burn than special 
fuels such as propylene oxide but that the high heat of combustion 
and availability of the petroleum hydrocarbons make it desirable 
to develop combustors for their use. 

The distribution of fuel in air is of importance in determining 
the performance of a combustor. It has been found that fuel 
sprays spread through a turbulent stream of air at a rate pro- 
portionate to the concentration gradient. With knowledge of 
this proportionality constant, it is possible to design an injector 
system to obtain the desired fuel distribution. 

The problems of flame stabilization and flame spreading are 
more complex, but it has been found that the stability of a flame 
held in the wake of a baffle is capable of some generalization, and 
it is felt that further investigation of the variables involved in 
these phenomena will lead, in time, to a sound design basis for 


tam-jet combustors 


INTRODUCTION 


es POTENTIALITIES OF THE RAM-JET aS a power 
plant for high-speed flight have led, in recent years, 
to its extensive investigation under the sponsorship 
of various Government agencies. The ram-jet seems 
best suited for supersonic flight, and its characteristic 
of requiring a high velocity to develop a practical 
amount of thrust and its simplicity of construction 
indicate that its most promising application may be 
in high-speed guided missiles where only a one-way 
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trip is required and the problems of take-off and land- 
ing are minimized. Its apparent lack of application to 
commercial usage, and the military importance of its 
probable applications prevent public discussion of 
much of the material now available on ram-jet per- 
formance and details of design; however, a brief dis- 
cussion of a few of the many technical problems in- 
volved in combustor development can be given. The 
bulk of the work discussed in this paper was done in 
cooperation with the Applied Physics Laboratory of 
The Johns Hopkins University under contract with 
the U.S. Navy Bureau of Ordnance. 

Because of the importance of frontal area in deter- 
mining the drag, and therefore the power requirement, 
at supersonic speeds it is necessary to keep the diameter 
of a ram-jet as small as possible for a given power out- 
put. This results in air velocities entering the com- 
bustion chamber sometimes greater than 300 ft. per 
sec. under normal operating conditions and higher than 
600 ft. per sec. under starting conditions for some launch- 
ing sequences. The combustion problems encountered in 
turbojet development due to high air velocity are there- 
fore accentuated, since a typical air velocity preceding 
the combustion chamber of a turbojet would be ap- 
proximately 50 ft. per sec. While it is permissible to 
allow 15 velocity heads of pressure drop in a turbojet 
combustor to obtain satisfactory combustion and mix- 
ing, this pressure drop would be disastrous at the higher 
velocities encountered in a ram-jet. The problem is 
eased somewhat in a ram-jet by the relative unimpor- 
tance of having a well mixed stream of hot gases leave 
the combustion chamber; however, the fact that in a 
ram-jet it is necessary to burn efficiently in a higher 
velocity air stream with less drag due to the combus- 
tion chamber has given rise to problems and designs 
differing considerably from those found in turbojet 
combustion chambers. 

The temperature and pressure in the combustion 
chamber depend on altitude and flight velocity and vary 
over a wide range depending on operating conditions 
Air temperatures between 100° and 1,200°F. might be 
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encountered and pressures range from several atmos- 
pheres down to a fraction of an atmosphere While 
the total range of conditions might not be encountered 
by any one ram-jet, the general problem of combustor 
development does cover such a range. 

Although the principle of the ram-jet was proposed 
some time ago in France, the successful application of 
the principle to flight was not made until 1945 when 
supersonic flight was maintained by a ram-jet developed 
by The Johns Hopkins University’s Applied Physics 
Laboratory and associated contractors under the spon 
sorship of the U.S. Navy Bureau of Ordnance. As in 
most combustion equipment, development was carried 
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out by trial and error, since the combination of the 
many important variables involved makes application 
of the principles of hydrodynamics and reaction kinet- 
ics difficult. Facilities for testing the larger sizes oj 
engines, however, are extremely expensive to construct 
and to operate, and engines of a number of sizes operat- 
ing under widely different conditions are required. [t 
is therefore important not only to develop engines by 
trial and error methods but to work out a set of basic 
design principles that will allow design of large units 
to give satisfactory operation with a minimum of modi 
fication and testing. 

Combustor development seems to a great extent, to 
be a matter of choosing a satisfactory fuel, mixing the 
fuel properly with the air, igniting it, and then allowing 
sufficient time for the combustion reaction to run its 
course. For purposes of the discussion that follows, 
the subject matter has been divided into the above 


categories. 


CHOICE OF FUEL 


Operational requirements determine. to a great ex- 
tent, the air temperature, pressure, and velocity, but 
some choice is allowed in fuel, which also is important 
in determining the design of the combustion chamber 
The ram-jet differs from the turbojet in this respect, 
since, while an airplane during its lifetime may consume 
a large quantity of fuel compared to its weight, the 
ram-jet, if used in a missile, makes only one flight and 
uses a relatively small quantity of fuel. So, while the 
problem of supplying sufficient quantities of fuel may 
control fuel specifications for piloted aircraft, missiles, 
because of the relatively small quantity of fuel used, 
could, if desirable, be fueled with more expensive fuels, 
which, perhaps, could only be produced in limited 
quantity. On the other hand, it is always desirable to 
use a cheap and plentiful fuel, so that any special fuel 
would have to show a distinct advantage. 

Aside from availability, a fuel would be judged by the 
1) heat of combustion per 


{ 


following characteristics: 
unit volume of fuel and per pound of air; (2) ease ol 
handling; and (3) ease of combustion. 

Ease of handling tends to rule out solid and gaseous 
fuels, since liquid fuels are much to be preferred by this 
criterion. Unstable fuels that tend to decompose or to 
become dangerous in storage would also tend to be ruled 
out. 

Of the stable liquid fuels presently available in rea 
sonably large quantities, petroleum hydrocarbons give 
the highest heat release per unit volume but are some 
what inferior to some oxygenated hydrocarbons in ease 
of combustion. This point is illustrated in Fig. 
(reference 1) where combustion efficiency is plotted as4 
function of equivalence ratio (stoichiometric air/fuel 
ratio divided by mixture air/fuel ratio) for a ram-jet 
type combustor burning pentane and propylene oxide 
It is seen that this particular combustor does not but! 
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COMBUSTION 
pentane so efficiently as propylene oxide. The propyl- 
ene oxide, on the other hand, has a gross heating value 
of 14,167 B.t.u. per Ib. of fuel as compared to 20,810 
B.t.u. per Ib. for pentane. Since propylene oxide is 
also more expensive and less plentiful than pentane, 
it is desirable to develop a burner capable of success- 
fully burning a hydrocarbon such as pentane. 

Several important points are brought out by con 
sidering measurements of laminar burning velocity and 
static inflammability limits. Fig. 2 shows the laminar 
burning velocity, determined by the Bunsen burner 
method, of a petroleum naphtha, kerosene, and pro- 
It is 
seen that the kerosene and the naphtha have nearly 


pylene oxide as a function of equivalence ratio. 


the same burning velocity but that propylene oxide 
will burn a great deal faster. Maximum burning ve- 
locity occurs slightly on the rich side of stoichiometric, 
indicating greatest ease of burning under rich condi- 
tions. It is also noted that, as temperature is reduced 
from 400° to 200°F., the maximum burning velocity is 
reduced approximately 24 per cent, so one would ex- 
pect more difficulty in developing a ram-jet to operate 
satisfactorily at low air temperatures than at high 
ones. It is also of interest to note that the highest 
flame velocity on Fig. 2 is 1.7 ft. per sec., which is small 
compared to the velocity of the air stream in which the 
flame is to be stabilized. 

Fig. 3 shows static inflammability limits as a function 
of pressure for several fuels. It will be noted that all of 
these fuels have a minimum pressure below which igni- 
tion cannot be obtained and that this minimum pres- 
sure occurs at mixtures somewhat richer than stoichio- 
metric. 

Iso-pentane has the narrowest range over which a 
flame will propagate and also the highest minimum ig- 
nition pressure. Propylene oxide can be burned at a 
lower pressure and over a wider range, and acetylene 
can be burned over a much wider range of conditions. 
It would be expected that, under the difficult burning 
conditions occurring in a ram-jet combustor, one would 
not quite reach the limits of operating indicated by the 
however, the static in- 
If the 
same improvement of minimum ignition pressure were 


static inflammability tests; 
flammability limits probably indicate trends. 


noted in a ram-jet combustor as in inflammability 
limits, then one might expect to be able to operate at 
an altitude approximately 35,000 ft. higher if acetylene 
were used rather than iso-pentane. This advantage 
would be difficult to realize in practice because of the 
difficulty in handling acetylene. 

The condition of the fuel (whether it is liquid or 
vapor) when it arrives at the point of ignition also has 
an effect on its combustion characteristics. It is 
therefore of interest to know whether a fuel, when 
injected into the combustion chamber, is capable of 
evaporating before reaching the point of ignition. It is 
possible to calculate the equilibrium fraction evapo- 
tated of a fuel for various conditions if the distillation 
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FIGURE 3 
STATIC INFLAMMABILITY LIMITS AS 
A FUNCTION OF PRESSURE 
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curve and the molecular weight of the various fractions 
distilled are known. Such information is available for 
petroleum fuels, and the results of such a calculation 
are shown in Fig. 4 for a stoichiometric mixture of a 
typical kerosene-type fuel and air under various condi- 
tions of pressure and temperature. The whole range 
of pressures and temperatures shown on Fig. 4 may be 
encountered in flights, and it can be seen that under 
some conditions of flight the kerosene is capable of 
completely evaporating, while under other conditions 
only a small fraction can evaporate. It would, of 
course, be possible to choose a fuel of such low vapor 
pressure that evaporation would be kept small or a fuel 
so volatile that substantially complete evaporation 
Such 
considerations are of importance in selecting a fuel 


would be possible under all operating conditions 


EQUILIBRIUM VAPORIZATION OF KEROSENE AS A FUNCTION 
OF PRESSURE AND TEMPERATURE 





STOICHIOMETRIC MIXTURE 


| FIGURES ABOVE CURVES ARE DIFFUSER 
EXIT TEMPERATURES 





PER CENT VAPORIZED 














.¢] 0.5 me 15 2.0 2.5 3.0 3.5 4.0 
PRESSURE - ATM. ABS 








710 JOURNAL OF THE 





FIGURE 5 


DIAGRAMMATIC SKETCH 
OF A RAM - JET 


FUEL — 


COMBUSTION CHAMBER 


py 5 a a 


DIFFUSER 
USE FLAME HOLDER 


Combustion, in general, proceeds more easily if the fuel 
is at least partially vaporized at the point of ignition; 
so it would be desirable to choose a fuel of volatility 
characteristics such that evaporation is possible under 
the conditions of flight. Handling characteristics, 
however, must also be considered, since fuels of ex- 
tremely low volatility tend to freeze under low-tempera- 
ture storage conditions and fuels of high vapor pressure 
tend to evaporate in the fuel tank or cause vapor lock 
in the pumping system under some conditions of high 
altitude operation. 

In view of the short time elapsing between the time 
a particle of fuel is injected and the time it reaches the 
point of ignition, it cannot be expected that equilibrium 
will be reached; however, such calculations seem to 
define the limits to which evaporation can go. The 
extent to which equilibrium is approached is consider- 
ably more difficult to determine and is a problem of 
current interest. 


FIGURE 6 
FUEL CONCENTRATION DOWNSTREAM FROM A SINGLE 
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Fig. 5is a diagrammatic sketch of a common type 
of ram-jet design. Fuel is injected upstream of 
the igniting device, which is followed by the com- 
bustion chamber in which burning and mixing are 
completed. These three sections will be discussed 
separately, since their problems and functions are 


different. 


MIXING OF FUEL AND AIR 


The purpose of the fuel injection system is to provide 
a well-atomized supply of fuel, properly mixed with the 
air, to the igniting section. The high-velocity stream 
of air occurring in a ram-jet is an advantage in obtain- 
ing good atomization since the impact of the air stream 
on the entering fuel is quite effective in breaking the 
fuel up into small particles. One could estimate from 
the work of Nukeyama and Tanesawa? that, with an 
air stream having a velocity of 300 ft. per sec. relative 
to the fuel jet, the weight median drop diameter would 
be approximately 40 microns. In a turbojet, such a 
high velocity is not readily available, and vortex-type 
mechanical atomizing nozzles are used. It would be 
estimated that with 100-lb.-per-sq.in. pressure drop 
through a typical turbojet vortex nozzle discharging 
into still air, a median drop diameter of 200 microns 
would be obtained. Since the vortex nozzle is satis- 
factory for many high-speed combustion uses, it is ap- 
parent that air atomization is promising for use in ram- 
jets. 

As illustrated in Figs. 1-3, the relative proportions 
of fuel and air are important in determining the ease of 
combustion. One would not expect flame propagation 
outside the inflammability limits, and greatest ease of 
burning would be expected for near stoichiometric 
mixtures. Since the average fuel/air ratio is deter- 
mined by operational requirements and may be con- 
siderably leaner than stoichiometric and since experi- 
ence has also shown that stability of operation is de- 
pendent on fuel distribution, it is important to be able 
to determine the distribution of fuel at the point of 
ignition. 

The distribution of fuel downstream from a single 
contrastream point of fuel injection mounted in a 6-in. 
diameter pipe is shown in Fig. 6. 

The data points were determined by measurement of 
the fuel/air ratio at various points along a pipe diameter 
at distances of 10 and 20 in. downstream from the point 
of injection. The solid lines are plotted from the fol- 
lowing equation: 


f = (Wyu/4nCW,X)e7 "4X 


where 
f = fuel/air ratio 
W, = weight rate of fuel flow 
W, = weight rate of air flow per unit of duct area 
u = air velocity 
r = distance from center of duct 





pl 
we 
mi 
eq 
air 
cel 


an 
ma 
ope 
fun 
vel 


tvpe 
m of 
com- 
g are 
‘ussed 
S are 


‘ovide 
th the 
tream 
otain- 
tream 
g the 
from 
th an 
lative 
would 
uch a 
-type 
Id be 
drop 
irging 
icrons 
satis- 
is ap- 
ram- 


rtions 
ase of 
ration 
ase of 
netric 
leter- 
' con- 
x peri- 
is de- 
> able 
nt of 


single 
| 6-in. 


ent of 
meter 
point 
e fol- 


irea 





COMBUSTION 


X = distance from point of injection to diameter 
along which samples were taken 
C = diffusion coefficient 


A diffusion coefficient of 0.705 sq.ft. per sec. was used, 
and it is seen that this equation describes the distribu- 
tion of fuel with sufficient accuracy to be adequate for 
many purposes. The above equation was derived by 
assuming that rate of diffusion was proportional to con- 
centration gradient and that injection point was a point 
A derivation of the 
equation is given in the Appendix. With knowledge of 
how this diffusion coefficient varies with other variables 
such as velocity, pressure, fraction vaporized, etc., it is 


source of the diffusing material. 


possible to calculate the fuel distribution for more com- 
plicated cases, such as when several injection points 
are used or when wall effects are important. 


FLAME STABILIZATION 


It has been found that under most ram-jet operating 
conditions a flame can be stabilized in the wake of a 
baffle mounted in the flowing mixture of fuel and air. 
This method of ignition is the basis of many ram-jet 
combustor designs and has been given considerable 
study. It is believed that the flame is stabilized be- 
cause of a zone of recirculation established in the wake 
of the baffle. By means of this recirculation, hot com- 
bustion products are fed back into the baffle where they 
serve aS a continuous source of ignition for the fresh 
mixture flowing past the edge of the baffle. Under 
high-velocity flow conditions with a baffle of circular 
cross section, material can flow upstream to the baffle 
from a point as far as 1.7 baffle diameters downstream 
from the baffle. Under adverse operating conditions, 
the flame is observed to pinch down at the downstream 
end of this recirculation zone and seems first to be- 
come unstable at this point. It is important to know the 
range of equivalence ratios over which a stable flame 
will be maintained in the wake of a baffle as a function 
of its size and shape and of the operating conditions. 
This was done in equipment in which some of the com- 
plicating factors that occur in a complete ram-jet 
were eliminated. Instead of an unevenly distributed 
mixture of fuel and air, as encountered in jet propulsion 
equipment, a homogeneous mixture of fuel vapors and 
air flowed past the baffle. The baffle was mounted 
centrally in the duct so that flow was symmetrical, and 
an effort was made to avoid large pressure fluctuations 
by using a short tailpipe. This equipment is diagram- 
matically illustrated in Fig. 7; in this case the baffle is 
of a long cylindrical shape mounted axially in the 6-in. 
pipe. 

The results of this study for 1 atmosphere pressure 
and an entering air temperature of 300°F. are sum- 
marized in Fig. 8. In this figure the range of stable 
operation in terms of equivalence ratio is plotted as a 
function of the baffle radius (R) divided by the mixture 
velocity (V) past the baffle. At high values of R/V 
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FIGURE 7 
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the curves become almost parallel, the lean limit ap- 
proaching a value of 0.63 and the rich limit approach- 
ing a value of 1.7 indicating that above an R/V of 20 X 
10~ sec. little gain in stability range is to be expected 
by using larger flameholders for a given mixture ve- 
locity. Substantial losses in stability range seem to 
occur, however, if baffle diameters and velocity are such 
that R/V is much below this value. It is interesting 
to note that the lean limit of 0.63 is very near the static 
lean inflammability limit of 0.58 for this fuel but that 
the rich limit of 1.7 is much leaner than the static in- 
flammability limit of approximately 3.4 generally re- 
ported for this type of fuel. 

Baffles of other shapes have been tested, and the re- 
sults can be summarized by stating that, for baffles of 
circular cross section, streamlining of the downstream 
face tends to reduce the stability of the flame but that 
other changes do not have a large effect by this cri- 
terion. It was also found that the results for baffles of 
other cross-sectional shapes can be placed on Fig. 8 
by using twice the mean hydraulic radius for the char- 
acterizing dimension. A similar study on smaller 
baffles was carried out at M.I.T.* The results were, 
in general, similar; however, they found that R/V 
gave the best correlation for their data. 

The effect of inlet air temperature is illustrated in 
Table 1, where it is seen that increasing the inlet air 
temperature tends to widen the limits, From meas- 
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FIGURE 9 
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urement of static inflammability limits one would have 
expected approximately the degree of widening ob- 


served, 


TABLE 1 
Effect of Temperature on Stability Limits 
1 Atmosphere; Air Velocity: 275 Ft. per Sec.; Baffle 
Diameter: 0.75 In. 


Pressure: 








Stability Limits 


Inlet Air Temperature, 
oF. Equivalence Ratio 


205 1 .27-0.77 
330 1.37-0.76 
400 1.37-0.73 

1.37-0.69 


500 


The effect of pressure is indicated by the data in 
Table 2. A larger effect at the low pressure was ob- 
served with smaller flameholders; however, with the 
larger flameholder, range of operation was quite in- 
sensitive to changes in pressure. 


TABLE 2 
Effect of Pressure on Stability Limits 
Inlet Air Temperature: 300°F.; Baffle Diameter: 2.87 In 
Range of Stable Burning 
Equivalence Ratio 


Velocity, 
Ft. per Sec. 


Pressure, 
Lbs. per Sq.In. Absolute 


2.9 +70 1.31-0.69 
15 400 1.68-0.60 


40 $00 1.68-0.61 


Using R/V as a criterion of stability it is possible to 
estimate the relative stability of various baffle configur 
ations. For example, if an extremely small baffle is 
placed in a duct through which a given mass flow of 
air flows, increasing its size a small amount will result 
in little change in the velocity past it compared to the 
increase in R, so a higher value of R/ |’ will be obtained. 
However, if the process is continued, the blockage of the 
baffle will become appreciable, the velocity will increase 
considerably, and the value of R/V will diminish. 
Fig. 9 shows R/V as a function of blockage ratio, and 
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NCES DECEMBER, 1949 
it is seen that the highest R J is obtained at a blockage 
ratio of 0.3. 

The above calculation is an example of the current 
attempt to develop design criteria for ram-jet com 
bustors; however, it is far from a complete study, since 
the data used were taken under idealized conditions, 
and other considerations of equal importance such as 
drag and combustion efficiency enter into flame sta- 


bilizer design. 


FLAME SPREADING 


The combustion chamber receives a mixture of air 
and fuel in which combustion has been started at a 
number of points, and it is necessary that sufficient 
space be allotted for the flame to spread through the 
fuel-air mixture. One would expect that the rate of 
combustion depends on the burning velocity normal to 
the flame front and on the area of the flame front. 

While the laminar burning velocity, determined by 
the Bunsen-burner technique, is measured under flow 
conditions quite different from those occurring in a 
rami-jet, data indicate that the laminar burning velocity 
will, under some conditions, serve as a measure of the 
relative amount of combustion of a fuel in a ram-jet. 
This was pointed out by Scurlock,* who compared, by 
time exposure photographs, the rate of spread of city 
gas-air and propane-air flames under simulated ram-jet 
operating conditions. He that the 
spread of the flame was greater for the city gas-air mix 


found rate ol 
ture, which had approximately twice the laminar burn 
ing velocity of the propane-air mixture, approximately 
in proportion to the ratio of laminar burning velocities. 
This effect is also illustrated by Figs. l and 2. In Fig. 
| it is shown that the combustion efficiency for a par- 
ticular ram-jet is highest for mixture ratios near stoichi 
ometric and that propylene oxide gives a higher com 
bustion efficiency than a paraffin hydrocarbon. Ex 
amination of Fig. 2 shows the same trend for the laminar 
burning velocities of these fuels, 

Since the thickness of the flame front under laminar 
burning conditions is of the order of 0.1 mm., the scale 
of turbulence existing in large equipment would be ex- 
pected to be considerably larger than this. In this case 
the area of the flame front is increased by the distortion 
of the surface due to turbulent flow. Photographic 
evidence presented by Scurlock® indicates that, under 
ram-jet conditions, the turbulence of the air upstream 
from the igniters has relatively little effect on the rate 
of flame spreading as compared to the mixing intro- 
duced by the flow past the igniters and the mixing 
caused by the acceleration of the burned gases. High- 
speed photographs of flame fronts show that the sur- 
face of the flame is extremely ragged. Detefmination 
of the area is therefore difficult, and the tendency has 
been to report a flame velocity based on time exposure 
photographs of the flame. Such a value would be high 
if used to calculate the rate at which the mixture is 
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burned, since a large amount of unburned fuel may be 
included within the boundaries indicated by such a 
photograph. Quantitative data on actual rates of 
burning are not available; however, this problem is of 
considerable theoretical and practical interest. 

Since the flame spreads through the unburned mix- 
ture from the surface of the flame held in the wake of the 
flame stabilizers, one would expect that the length of 
combustion chamber would then depend on the distance 
between flame stabilizers. In order to keep this dis 
tance to a minimum, one would want to use a large 
number of flame stabilizers, but to keep the blockage 
of the stabilizer low it would be necessary to use small 
stabilizers. In the preceding section, however, it was 
pointed out that as the size of a baffle is reduced the 
stability range is also reduced; so it is necessary to 
compromise between the shorter combustion chamber 
possible with a large number of small baffles and the loss 
of igniter stability resulting from reducing the baffle 
size. In order to arrive at the best compromise for the 
great variety of possible ram-jet motors, it can be seen 
that accurate data on baffle stability and on combus 
tion rates are necessary. 

One of the outstanding problems encountered when 
the combustion chamber is added to the flame stabiliza- 
tion section is the rough and explosive burning that 
sometimes occurs. Knowing the fuel distribution and 
the stability range of the baffles, one could attempt to 
predict the range of air/fuel ratios over which stable 
operation can be maintained; however, rough and un 
stable burning sometimes occurs within this predicted 
range This phenomenon is apparently a consequence 
of spreading the flame throughout the high-velocity 
gas stream and might be expected from the following 
reasoning. Since the pressure at the inlet to the com 
bustion chamber can, under some circumstances, be 
twice that at the exit, because of the heating of the gas 
stream, a change in the amount of combustion will 
change the pressure, and this change in pressure can 
change the flow through the diffuser, the rate of fuel 
flow, and other variables of importance to the combus 
tor. With such a situation many opportunities for 
resonance phenomena and instability exist, depending 
on the exact relation between these variables, and, while 
rules of thumb that allow one to avoid this situation are 
being developed through experience, this is another 
phenomenon that needs study and generalization 


CONCLUSIONS 


While the data and discussion given in this paper 
have perhaps raised more questions than they have 
given answers, it is felt that indications exist that gen 
eralization of many of the processes involved in ram-jet 
combustion is possible. Development of a firm basis 
for the various phases of combustor design will require 
a large amount of further research, which, it is felt, is 
justified in view of the great expense of testing large 
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ram-jets aud the variety of conditions under which it 
appears ram-jets will operate. 
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APPENDIX | 


Derivation of Equation for Fuel Concentration 
Downstream from a Point Source 


If it is assumed that the rate of spreading of fuel in 
a ram-jet downstream from a point of fuel injection is 
proportional to the concentration gradient, it is possible 
to calculate the shape of the fuel distribution curve as 
follows: 


O1'/Ofr CA (Of /Or) (1) 
where 
|’ = the volume of diffusing material 
{ = time 
C = diffusion coefficient 
f = fuel/air ratio 
A = area normal to the direction r 


Chis equation can be applied to the case where fuel 
is injected at a point making the following assump- 
(1) steady flow in the system; (2) a flat velocity 
distribution; (3) the amount of mixing in the direction 


tions: 


of flow (xv) negligible in comparison with the mixing in a 
radial direction (7); (4) the mixing coefficient (C) con- 
stant at ail points in the svstem; and (5) the duct is 


extremely large in diameter. 


\pplying a material balance to this system, Eq. (1) 


becomes: 
Of/ox = (C/u)[(1 (Of Or) + (0°f/Or?)| (2) 
where “ is the stream velocity. 


\ solution of the equation which will also describe 


the boundary conditions is 


f = (l/ax)e - (3) 


where u and # are constants 


By differentiating Eq. (3) and substituting in Eq. (2), 
one finds that 


) u/AC (4) 
Substituting Eq. (4) into Eq. (3), we have 


f = (l/ax)e~“"" *“* (5) 


(Continued on page 760) 








Basic Studies on Flame Stabilization 


GLENN C. WILLIAMS* 


Massachusetts Institute of 


ABSTRACT 


This is a report of initial studies to determine the mechanism of 
anchoring flame in high-velocity air-fuel streams. Simple bluff 
shapes were used as flame stabilizers in a constant area 1- X 3-in. 
rectangular duct carrying homogeneous mixtures of air and 
gaseous fuel. The limiting gas velocity above which flames blew 
off the stabilizers were studied as it was affected by: 

(1) The approach gas stream variables. (a) Fuel-air ratio (over 
the whole range of inflammability). (b) Turbulence (intensities 
from 0.3 to 8 per cent of the approach stream velocity). (c) Fuel 
type (commercial propane and city gas). 

(2) Stabilizer variables. (a) Size (Characteristic dimension 
from 0.02 to 0.5 in.) (b) Shape (simple, faired, and shortened 
rods, flat plates, gutters, spheres). (c) Addition or extraction of 
heat from stabilizer. 

On the basis of experimental evidence, the following con- 
clusions were reached: 

(1) Flame stabilization behind bluff obstacles depends on the 
presence of a recirculating eddy region that is filled with hot 
combustion products and furnishes continuous ignition to the 
approaching unburned gases. 

(2) Over the range of stabilizer sizes studied, the blowout 
velocity under steady flow conditions was found to be affected 
primarily by the properties of the approach gas stream and the 
stabilizer thickness. 

(a) For span-thickness ratios greater than 2, the blowout 
velocity at a given air-fuel ratio in the absence of approach 
stream turbulence was independent of the stabilizer shape and 
proportional to the thickness to the 0.45 power. 

(b) For span-thickness ratios less than 2, the corresponding 
blowout velocity was roughly proportional to the first power of 
the thickness. 

(c) The stability limits for a stabilizer of given characteristic 
dimension are unaffected by variation (simulated) of ratio of 
chamber width to stabilizer thickness over the range 10 to 80. 

(3) Increasing turbulence in the approach gas stream decreased 
both the velocity and air-fuel ratio ranges over which flame could 
be maintained. The effect of imposing turbulence of a given scale 
and intensity diminishes with increasing ratio of stabilizer thick- 
ness to turbulence scale. 

(4) Heating the stabilizer other than by the flame markedly 
increases the stability limits, and externally cooling the stabilizer 
markedly decreases the stability limits. 

(5) The maximum blowout velocities occur with slightly leaner 
than stoichiometric mixtures of city gas (26 per cent He) and 
richer than stoichiometric for propane. This is attributed to the 
different ratios of molecular diffusivity to that of air for the two 
fuels. 

(6) Flame blowoff usually proceeds stepwise, flame propaga 
tion throughout the main stream ceasing before the eddy zone be- 
hind the stabilizer is quenched 


NOMENCLATURE 


dimensionless constants 
mass specific heat at constant pressure 


a,c = 
Cp = 
Presented at the fourth National Flight Propulsion Meeting, 


I.A.S., Cleveland, March 18, 1949. 
* Professor of Chemical Engineering 
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Technology 


= mass ratio of fuel to air 

stabilizer thickness 

= absolute temperature 

velocity 

thickness of preignition zone 

= kinematic viscosity 

subscript indicating entrance conditions 
= subscript indicating ignition conditions 


(I) INTRODUCTION 


O' THE PERFORMANCE PARAMETERS used _ to 
evaluate combustion chamber performance for 
aircraft propulsion application, one of the most critical 
is the limitation of operating ranges of air-fuel ratio and 
flow velocity due to inability to maintain stable com- 
bustion. In order to attack the problem of designing 
combustion chambers having desirably wide ranges of 
operability, it is in order to attempt to discover the 
laws governing the mechanism of flame stabilization or 
at least to determine empirically the variables affecting 
flame stability in high-output systems, as well as quan 
titative information on the effects of variation of these 
variables. It is well known®* ° that, at the gas velocity 
ranges of interest in the operation of Bunsen-type 
flames, the criterion for flame stabilization on a surface 
parallel to the gas flow calls for the existence of a par- 
ticular gas velocity profile in a plane normal to the 
surface. This profile must include, at a distance no 
less than a certain minimum value (the quenching dis 
tance) from the stabilizing surface, a gas velocity no 
greater than laminar flame velocity. Since laminar 
flame velocities for most fuel-air mixtures of general in- 
terest are less than a few feet per second and since 
laminar boundary layers on surfaces in high-velocity 
streams are usually much thinner than the minimum 
allowable quenching distance, flame-holding on flat 
surfaces is obviously unsuited to high-output systems 
that may operate with inlet gas velocities as high as 
500 ft. per sec. Early experiments’~’ with flames 
stabilized downstream 
that flames could be stabilized on these wires at higher 
In contrast to ex- 


of blunt-ended wires showed 
gas velocities than on flat surfaces. 
periments on normal Bunsen flames, no correlation ap- 
pears to exit between the gas velocity at which flame 
blew off the blunt ends of the wires and the velocity 
gradient at the wire surface. The most successful 
steady-flow propulsive combustion chambers developed 
to date have exploited the superior flame-holding 
abilities of bluff shapes. Attempts to isolate separate 
effects of the various variables presumed to influence 
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burner performance by studies on complete combustion 
chambers have been unsuccessful. 

The purpose of this investigation was to study flame 
stabilization on bluff shapes in experiments designed to 
allow independent study of several variables while 
eliminating or holding constant those variables not 
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under study. The variables studied were air-fuel ratio, 
approach gas velocity, turbulence in the approach gas 
stream, flame stabilizer size and shape, and fuel type 
(Cambridge city gas and commercial propane). 


(Il) APPARATUS AND PROCEDURE 


The experimental apparatus, Fig. 1, consisted of a 
constant cross section, 3 by 1 in., low turbulence duct, 
17 in. long, through which homogeneous mixtures of 
gaseous fuel in air could be flowed steadily with en- 
trance pressure and temperature of about one atmos- 
phere and 300°K. and at entrance velocities from 20 to 
350 ft. per sec. 

As shown in Fig. 1, the mixture flows through a 
calming section—consisting of a diffuser, five 200-mesh 
screens, and a 25 to | contraction-ratio nozzle—and 
enters the combustion chamber. The turbulence in- 
tensity at entry is less than 0.5 per cent of the mean 
velocity, and the velocity profile is flat. Screens of 
different mesh and wire diameter were placed at the 
combustion chamber entrance to raise the turbulence 
level in the chamber. The rectangular combustion 
chamber had two 3-in. wide opposing walls of Vycor 
glass to allow observation and photography of the flame. 
Flame stabilizers of simple shape and of various sizes 
were supported with axes normal to the direction of 
flow between the opposing glass sidewalls. The stabil- 
ized flame when viewed along the axis of a stabilizer 
gave an appearance of two-dimensionality. 

The parameters chosen for study were varied as 
follows: air-fuel ratio, lean limit to rich limit blowout; 
combustion chamber entrance velocity, 20-350 ft. per 
sec.; entering stream turbulence intensity, 0.4—S.0 per 
cent; turbulence scale, 0.01—0.08 in.; flame stabilizer 
thickness, 0.016-0.5 in.. span, 0.1-1.0 in.; flame 
stabilizer shapes, single cylindrical rods, shortened rods 
with rounded ends, rods with wedge-shaped upstream 
fairings, spheres, 30° closed-edge and open-edge gutters, 
flat plates, multiple rods. 

Characteristics of the fuel-gases used were: 


Cambridge City Gas Commercial Propane 


CO:.. 7.3% C3H¢ 2.4 
Iluminants = 9.8% C;Hes 83.8 
Os» 0.6% CoH. 13.3 
CO : 21.8% CyHyw 0.5 
CH, 14.7% 
H, 26.8% 
Ne. . ee 19.0; 
Pounds air per pound gas required to make stoichiometric 


ih kata eee 15.6 


mixtures 6.51....... 


Runs were made to obtain blowout limits and photo- 
graphs of the flame at various air-fuel ratios and com- 
bustion chamber entrance velocities, with different 
sizes, shapes, and locations of the flame stabilizer, and 
for different turbulence levels. 

Prior to making runs with flame, the gas flow through 
the chamber was examined in some detail, with tur- 
bulence generated by one of three screens of mesh 3, 8, 
and 55 per in.; wire diameter, 0.063, 0.025, and 0.007 
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in.; and fractional blocked area, 0.542, 0.360, and 0.629, 
Figs. 2 
and 3 show the results of calculations and measurements 
of turbulence. The were 
Dryden’s equation? for decay of turbulence down 


respectively, as well as with no screen present. 


calculations made _ using 
stream from screens, assuming in the case of the no 
screen runs that the turbulence was that existing 6 in. 
from a 200-mesh screen (in the calming section), re 
duced twenty-five-fold in intensity by contraction but 
unchanged in scale. Experimental measurements of 
turbulence were made with a two-probe hot-wire 
anemometer designed for use up to frequencies of 20,000 
cycles per sec. The measured per cent turbulence falls 
below the calculated values for all three screens except 
near the exit end of the chamber, where the measured 
values begin to increase because of diffusion out into 
the main stream of turbulence generated in the bound 
ary layer. 

With no screen at the chamber entrance the measured 
values of turbulence are considerably above the calcu 
lated values (0.4 per cent measured, 0.02 per cent cal 
culated). That this discrepancy is probably due to 
pulsations from the air compressor is shown by com- 
parison of the 0.07 per cent turbuleace measured with 
the air supplied to the chamber by blowing down a large 
surge tank, with a reading 0.3 per cent obtained under 
the same conditions except with the compressor on. 
The measured values of scale are in every case larger 
than the calculated values. 
and calculated values of scale and intensity in Figs. 2 
and 3 suggests that the hot wire was not sufficiently 


Comparison of measured 


compensated for time lag in these tests. Other possible 
reasons for disagreement are the presence of small oil 
droplets and the vibratory effects of the air compressor. 

The critical velocity for boundary-layer turbulence 
was determined with the turbulence probe at three 
locations along the chamber length. Corresponding 
Reynolds Numbers based on the distance from the end 
of the converging section of the nozzle were about 10,° 
independent of location. At all velocities greater than 
about 20 ft. per sec., the boundary layer in the major 
portion of the chamber is turbulent with a measured in- 
tensity of about 3 to 5 per cent. No systematic change 
in boundary-layer thickness or ratio of local to mean 


velocity Vo was observed with changes in Vo. 


(III) THEORETICAL CONSIDERATIONS 


As the cold 
bluff shape, separate from it, and then contact the hot 
gases in the eddy region, heat and chain carriers are 


approach stream gases pass around a 


transferred from the eddy zone to the portion of the 
This 
cold gas soon reaches its ignition temperature and burns. 
During the time that this recently ignited gas parcel 
flows past the rest of the eddy region, heat and chain 
carriers flow back into the eddy region and replenish the 
previous depletion. The recirculation of the hot gases 


cold unburned stream nearest the eddy region. 
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in the eddy region thus furnishes a continuous source of 
ignition for the unburned gases. Assuming thermally 
controlled ignition, the required heat supply rate for 
ignition of the approach stream increases with Io, with 
the difference (7; 70) between the ignition tempera- 
ture and the approach stream temperature, the thick- 
ness of the preignition zone 6, and the heat capacity per 
unit volume pCp. The rate of heat supply necessary to 
assure continuous ignition is then proportional to 


Vol6(pCp)(T; — To)] 1) 

If 6, is assumed proportional to the ratio of thermal 

diffusivity k/pCp to laminar flame speed S,, relation | 
becomes 


Vok/Szt (T: — To) (2) 
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From heat-transfer considerations the rate of heat 
supply from the eddy region (size assumed proportional 
ton, the stabilizer thickness) to the unburned gas prior 
to ignition is proportional to the product of the transfer 
area, the coefficient, and the temperature difference 
between burned and unburned gases, 


[(nVo/v)(k n)|n(T, To) (3) 


The bracketed term in relation (3) allows for the factors 
affecting the heat-transfer coellicient (v is the kinematic 
viscosity of the gases), and the other two terms allow 
for area and temperature difference. For blowout to 
occur, the heat supply rate must be jnfinitesimally less 
than the heat required. Equating relations (2) and (3) 
and rearranging so that all quantities thought to be a 
function of air-fuel ratio appear on the right-hand side, 
the following relation results: 


vm! [ 3<e ay be (*) (4) 
te re o\, | 


where the new exponent a replaces ¢ (1 — c) and f is the 
mass fuel-air ratio. 

It is interesting to note the only effect of pressure on 
the blowout limits as predicted by Eq. (4) is that on 
density through its effect on the kinematic viscosity p, 
which varies inversely with the pressure. 
blowout velocity should be reduced as pressure is de- 


creased. 


(IV) ExXxPpreRIMENTAL RESULTS 


Visual observation of the flames leads to the con- 
clusion that the eddy region just downstream of the 
stabilizer consists of completely burned gases. 

Examination of flame photographs similar to those of 
Fig. + shows no effect of air-fuel ratio or of approach- 
stream velocity on the location with respect to the 


stabilizer of the first appearance of reaction. 


T—~ 48.9 ¥ ——~ 49.9 
i/t —~ 12,60 iff -- 10.42 


Vi 54.9 









v——— 51.6 


v ---- 110.6 


V ---- 216.3 
Aff — 8.00 i/¢ 4.66 i/f -- 7.82 





T—— 56.6 
/f— 5.54 Ife — 4.78 


v —-- 56.8 
i/t — 4.09 





Fic. 4. Typical flames on rod flameholders 


Thus, the 
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Inability to hold flame in a high-velocity stream on 
an unheated surface that does not provide an eddy re- 
gion was demonstrated by failure of all attempts to 
stabilize flame on a '/2-in. long flat plate lying parallel 
to the flow. 


Mechanism of Blowout 


It is postulated that flame blows out when the eddy 
region regains less heat from the newly ignited gases 
than is required to ignite these gases. At constant 
approach velocity, as the mixture ratio is changed to re 
duce the flame temperature, the average temperature of 
the eddy region drops, finally reaching so low a tem- 
perature that continuous ignition is no longer maintained 
and the flame dies. Flame does not immediately cease, 
however, but continues for a few revolutions of the 
eddy region until it cools sufficiently to be impotent. 

Examination of high-speed movies of blowout shows 
that the last region where flame exists is near the rod, 
indicating that the flame extinguishes by first ceasing to 
propagate into the main stream and then shortening 
until the eddy region temperature becomes too low for 
ignition. 

Considerable quenching occurs when the burned gases 
just to the rear of the eddy region are rapidly accelerated 
as they move farther downstream in a constant area 
duct. This quenching action must account for the 
failure of the residual flame to propagate into the main 
stream. It appears that the propagation of the flame 
above the stabilizer into the unburned gases (which 
fails to take place in the case of the residual flame) and 
the stabilization mechanism itself are independent of 
one another. Residual flames appeared in only a few 
instances—-on the lean side of stoichiometric when city 
gas was used as the fuel and on the rich side with pro 
pane. With the fuels and stabilizers employed in these 
tests it appears that failure of the flame to stabilize, 
rather than its failure to propagate, is in most instances 
the controlling phenomenon. Heating the stabilizer 
should tend to make propagation controlling; this was 
verified. 

Blowout limits of city-gas flames stabilized on nine 
rods varying thirty-fold in diameter from 0.016 to 
0.498 in. were obtained both with and without a screen 
at the chamber entrance. Fig. 5 presents the results, 
which are seen to be consistent with Eq. (4) except in 
the range of extremely lean mixtures where the effect of 
rod size reverses. This effect is thought to be due to the 
increased cooling of the eddy region by the large stabi- 
lizers; experiments with water-cooled rods support this 
view. 

In contrast with the small effect of approach-stream 
turbulence on propagation reported elsewhere,* the 
reduction in stability is pronounced (Fig. 6) excepton the 
two largest rods. It is believed that the effect of tur- 
bulence on stability limits is primarily due to the in- 
tensity and size of the gusts characterizing the variation 
from the mean value of the local velocity at the stabi- 
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(FUEL: CAMBRIDGE CITY GAS — ENTRANCE PRESSURE- | ATM, ENTRANCE TEMP 300 340°k) 
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lizer. Turbulence of a given scale and intensity should (Vg,,)°-*® as ordinate are shown in Fig. 7. The correla- 
then have less influence on the stability limits of large tion is satisfactory, but, as expected, the points scatter 
stabilizers. The data confirm this. somewhat on the lean side at low velocity because of the 

To test Eq. (4), a best value of the exponent a, 0.45, quenching action of the stabilizer itself, which is not 
was found by plotting Vs, vs. m to logarithmic coor- taken into account in the heat balance on which Eq. (4) 
dinates. Blowout limit data for rods replotted with _ is based. 
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0.032-in. tube stabilizer. (Fuel: city gas.) 


The effects of heating and cooling the stabilizer are 
shown graphically in Fig. 8 for city gas-air mixtures. 
Blowout velocities obtained with a 0.038-in. tube main- 
tained at 1,200°K. match those for a 0.126-in. unheated 
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rod at an air-fuel ratio of 9.5. As the mixture ratio 
approaches stoichiometric, the unheated rod rises in 
temperature, and the effect of the artificial heating be 
comes less. As expected, blowout limits for the un 
heated tube were wider than for an unheated rod of 
equal diameter. Similar results were obtained with 
propane-air mixtures. ! 


Stabilizer Shape 


Closed-edge-gutter blowout-limit data for city gas-ait 
mixtures were obtained with five gutters ranging in 
width from 0.025 to 0.5 in. Comparison of the stability 
limit curves for closed-edge gutters with similar curves 
for rods shows that, although they both have the same 
general shape and appearance, lean-side limits are 
leaner for the gutters. This is due, it is presumed, to the 
relatively smaller cooling action of the gutters on the 
eddy region. 

It has been already been shown that the blowout 
limits for flames stabilized on rods and closed-edge 
gutters were approximately the same. Blowout limit 
data were obtained with a rod, a closed-edged gutter, 
open-edge gutters with a ratio of total area to open area 
of 100 and 25, and a flat plate, all having approximately 
the same characteristic dimension of 0.5 in. That a 
single curve be drawn through all these data without 
great dispersion is shown in Fig. 9. On the basis of the 
blowout limit data it is concluded that the length of the 
eddy region behind gutters, rods, and plates of the same 
characteristic dimension is approximately the same. 


Multiple Rod Stabilizers 


Blowout limit data obtained with 1, 2, 4, and 8, 0.038 
in. rod stabilizers placed in a single plane normal to the 
flow in the chamber are plotted in Fig. 10. 
curve fits all except the S-rod data. The effect of 


A single 


changing the ratio of chamber width to stabilizer width 
over a simulated eightfold ratio is shown to be sub 
stantially negligible. Blowout as determined by 
visual observation occurred simultaneously from all 
rods. 

Under similar conditions the stability of propane-air 
flames on rods was found to be,considerably less than 
that of city gas-air flames, a result expected from the 
lower value of laminar flame velocity. The marked re 
duction in stability caused by the introduction of 
approach-stream turbulence was again noted with pro 


pane as fuel. 


Effect of Stabilizer Span-Thickness Ratio 


Longwell’ obtained blowout data with naphtha-air 
mixtures flowing in the range of gas velocities from 300 
to 1,000 ft. per sec., with flames stabilized on axially 
symmetrical flameholders axially located in a cylindrical 
duct 6 in. in diameter. His blowout data appear to 
correlate reasonably well either as a function of stabil 
izer characteristic dimension to the 0.45 power or to the 
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lated variation of chamber width 


first power. Because of the existence of Longwell’s data 
for larger size flameholders than those studied in this 
investigation and because doubts were cast upon the 
approach to two-dimensionality of the present experi 
ments, a program of investigation! of this phase was 
initiated in which tests were made of stability limits in 
propane-air mixtures of rods progressively shortened 
from the full l-in. depth of the test chamber to a span 
thickness ratio of unity. Each of the rods shorter that 

in. im span was fitted with hemispherical ends and 
supported by a wire of small diameter passing through 
the axis of the rod. As shown by Fig. 11, the effect on 
blowout velocity of shortening the span of 0.1-in 
diameter rods is negligible until the ratio of span to 
thickness is less than 2. Similar results were obtained 
with 0.18S-in. diameter rods. Data obtained with 
spheres are shown in Fig. 12 to give reasonable cor 
relation with the first power of the sphere diametet 
as first proposed by Longwell. It should be noted that 
blowout limits for rods proportional to the 0.45 power ol 
rod diameter and blowout limits for spheres propor 
tional to the first power of sphere diameter call for : 
sphere of a large size to produce more stable flames than 
If it is visualized that the 
a rod having a length 


a rod of the same diameter. 
superior stability limits of 
diameter ratio greater than 2 are due to the ability of the 


rod to sustain flame when a quenching gust of size equal 


to rod diameter passes across the eddy region, whereas « 
similar sized quenching gust passing across the edd) 
region of a sphere would leave no fire to restart the edd) 


region, it does not appear that the relative exponents ou 
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propane.) 


size for spheres and rods are valid outside the range of 
experimental variables covered here. 


Stabilizer Boundary-Layer Thickness 


The stabilizer used in initial boundary-layer studies® 
consisted of a rod mounted with its axis normal to the 
direction of gas flow and with a wedge pointing up- 
stream attached to its upstream surface. Wedge fair- 
ings of various lengths were used to control the bound- 
ary layer on the stabilizer at the point of entry of the 
cold approach gases into the preignition and combus- 


Effect of span-thickness ratio on stability limits. (Fuel: 


tion zones. The length of the wedge fairing was in- 
creased progressively to increase boundary-layer thick- 
ness, and long fairings were used to deliver a turbulent 
boundary layer to the combustion zone. The faired 
section was externally cooled and insulated from the hot 
portion of the stabilizer in an attempt to find only 
those effects due to boundary-layer variation. Un- 
fortunately, the scatter of the blowout data taken in 
this series of tests was so great as to mask any effects of 
boundary-layer variation. It can be concluded that, if 
the effects of changes in wedge fairing length were con- 
fined to changes in boundary-layer thickness and tur- 
bulence alone, the effect of change in boundary layer is 
not nearly so pronounced as that of changing stabilizer 
temperature. It can be argued that changing the length 
of wedge fairing not only changed the boundary-layer 
thickness and turbulence at entry to the eddy zone but 
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also changed the approach angle of the boundary layer. 
That the effect of such a change in approach angle is 
probably small is militated for in the negligible effect on 
stability shown by changes from gutter to rod to flat 
plate in Fig. 9. 


Molecular Diffusion 


Most of the observations of the characteristics of 
flame stabilization on bluff shapes can apparently be 
explained on the basis of a thermal ignition theory and 
the gross aerodynamics of the system. Failure to 
correlate blowout velocity data for different fuels in 
terms of the laminar flame velocity, however, and some 
visual observation on differences in the character of 
flames of city gas and of propane lead to the tentative 
conclusion that molecular diffusion plays a role of some 
importance in determining the air-fuel ratio range over 
which stable flames are obtained. The appearance of 
residual flames only on the rich side of stoichiometric for 
propane-air mixtures and only on the lean side of 
stoichiometric for city gas-air mixtures argues in favor of 
this hypothesis. Since the molecular diffusivity of pro- 
pane is less than that of oxygen and since those of hy- 
drogen and methane (the major combustible constitu- 
ents of the city gas) are greater than that of oxygen, it 
can be argued that in the case of propane the mixture in 
the eddy region is always leaner in fuel than that of the 
main stream and in the case of city gas the eddy region 
is always richer in fuel than the main stream. This 
thesis is supported by the fact that the peaks of blowout 
velocity curves for propane-air mixtures always appear 
on the fuel-rich side and those of city gas-air mixtures 
always slightly on the fuel-lean side, plus the fact that 
the color of the eddy region of propane flames was al- 
ways characteristic of a leaner mixture than that of the 
burning main stream. It would be expected, in the ab- 
sence of concentration gradients of fuel-air ratio, that 
maximum blowout velocity would occur for stoichio- 
metric mixtures, since maximum flame velocity and 
maximum heat generation per unit volume are ex- 
pected near the stoichiometric ratio. 


(V) CONCLUSIONS 


(1) Flame stabilization in high velocity, premixed gas- 
eous fuel-air streams is not possible on unheated surfaces 
that do not produce recirculating eddy zones. Flame 
stabilization in the wake of an eddy-producing stabilizer 
depends on the presence of an eddy region filled with 
combustion products sufficiently hot to furnish contin- 
uous ignition to the approaching unburned gases. 

(2) Over a range of stabilizer sizes from 0.016 to 0.50 
in. in thickness and from 0.1 to 1.0 in. in span, the blow- 
out velocity for a given fuel in air was affected primarily 
by the gas stream variables of air-fuel ratio and tur- 
bulence and by the stabilizer thickness: 

(a) For span-thickness ratios greater than 2, the 
blowout velocity at a given air-fuel ratio and stream 
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turbulence was independent of stabilizer shape- 
cylindrical rods, rods with wedge-shaped upstream 
fairings, hemispherically ended cylindrical rods, 
open-edged gutters, flat plates. The blowout velocity 
was proportional to the 0.45 power of the stabilizer 
thickness. 

(b) For span-thickness ratios less than 2, the 
corresponding blowout velocity was roughly pro- 
portional to the first power of the thickness. 

(c) The stability limits for a stabilizer of given 
thickness are unaffected by simulated variation from 
10 to 80 of the ratio of stabilizer thickness to chamber 
width. 

(d) Increasing turbulence in the approach gas 
stream decreases the range of approach velocity and 
air-fuel ratio over which flame can be maintained 
The effect of approach stream turbulence of a given 
scale and intensity diminishes with increasing ratio of 
stabilizer thickness to turbulence scale. 

(3) Heat addition externally to a stabilizer holding 
flame in a gaseous combustible mixture markedly widens 
the stability limits. Removal of heat from such a 
stabilizer narrows the stability limits. 

(4) Flame blowoff from a bluff stabilizer usually pro 
ceeds stepwise. Flame propagation downstream from 
the eddy zone usually ceases before the eddy zone itself 


is quenched. 


(5) Shifting of the peak velocity of stable burning 
away from the stoichiometric fuel-air ratio normally ex- 
pected to give maximum blowout velocity is attributed 
to the mechanism of molecular diffusion operating to 
change the fuel-air ratio in the eddy zone from that in 
the main gas stream. 

(6) The effect on stability limits of changing fuel type 
is in the direction of the accompanying change in 


laminar flame speed. 


REFERENCES 


1 Baddour, R. F., and Carr, L. D., S. M. Thesis in Chemical En- 
gineering, Massachusetts Institute of Technology, 1949. 

2 Dryden, H. L., Turbulence and Diffusion, Ind. Eng. Chem., 
Vol. 31, pp. 416-425, April, 1939. 

3 von Elbe, G., and Mentser, M., Further Studies of the Struc- 
ture and Stability of Burner Flames, J. Chem. Phys., Vol. 13, pp 
89-100, February, 1945. 

4 Goss, W. H., and Cook, Emory, The Ramjet as a Supersoni 
Propulsion Plant, Trans. S.A.E., Vol. 2, No. 4, pp. 642-657, 
October, 1948. 

5 Lewis, B., and von Elbe, G., Stability and Structure of Burner 
Flames, J. Chem. Phys., Vol. 11, pp. 75-97, February, 1943. 

6 Longwell, J. P. (see reference 4), Standard Oil Development 
Company. , 

7 Mache, H., Die Physik der Verbrennungserscheinungen, \eit 
und Comp., pp. 25-54, Leipzig, 1918. 

8 Scurlock, A. C., Se.D. Thesis in Chemical Engineering, Massa 
chusetts Institute of Technology, 1948; abstracted in Flame 
Stabilization and Propagation in High-Velocity Gas Streams, 
Meteor Report No. 19, Fuels Research Laboratory, Massa- 
chusetts Institute of Technology, July, 1948. 

® Kellogg, W. W., S.M. Thesis in Aeronautical Engineering 
Massachusetts Institute of Technology, 1949. 





shape- 
stream 

rods, 
le city 
bilizer 


2, the 
pre )- 


given 
1 from 
unber 


h gas 
Vv and 
ained 
given 
itio of 


ding 
videns 
ich a 


"4 pro- 
from 
itself 


rning 
ly ex- 
duted 
ng to 
lat in 


type 


ye ill 


al En- 
‘hem., 


OSfric- 


3, pp 


rsonu 


57, 


urner 
3. 

nent 
, Vert 
[ussa 
Flame 
cams, 


[assa- 


ring, 





Equations for the Approximate Solution of 
Dynamic Problems for Stable Linear Systems 


KURT H. HOHENEMSER* 
McDonnell Aircraft Corporation 


SUMMARY 


In generalization of the ‘‘square wave"’ formulas, presented re- 
cently in this magazine,! a number of equations are derived which 
ire applicable in the approximate analysis of complicated linear 
systems or in the analysis of dynamic test data of a linear system. 
A summary of the obtained equations is given in the Appendix 


(1) INTRODUCTION 


I THE ANALYSIS AND SYNTHESIS Of linear systems, for 
example of aircraft automatic controls, it is often 
desirable to use approximate numerical or graphical 
methods for the calculation of a transient if the fre- 
quency response of the system is given or for the calcu- 
lation of the frequency response function if a transient 
is given. Approximate methods are indicated if either 
the exact analytical solution is too involved (for ex- 
ample, in systems where the characteristic equation 
has more than four roots) or if the system is given by 
measurements of the frequency response or of a transient 
rather than by its physical constants. In reference 1 
it was shown how well the frequency response function, 
calculated by an approximate method from the meas- 
ured airplane motion after an elevator step displace- 
ment, fitted the frequency response function as meas- 
ured in flight. The advantage of substituting the la- 
borious measurements of the airplane response to har- 
monic control displacements over the entire necessary 
frequency range by a single measurement of a transient, 
following a step displacement of the control, is obvious. 

It is often of importance to know the relations be- 
tween transient and frequency response in case of a 
gradual transition from one position of the control to 
another position with finite rate of change of the stick 
position. This procedure may be called oblique step 
displacement of the control, contrary to the rectangular 
step displacement. Finally, an alternative method 
of exciting a transient of an airplane is the use of an 
impulse control deflection instead of a step deflection. 
The impulse method consists of deflecting the control 
for a short time and going back to the original position. 

In the following, a number of approximate formulas 
is developed appropriate for numerical or graphical 
determination of the frequency response from a given 
transient and of a transient from a given frequency 
response covering the two cases of oblique step input and 
of impulse input. The case of rectangular step input, 
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treated in reference 1, is included in the more general 


equations for oblique step input. 


(2) APPROXIMATE CALCULATION OF A TRANSIENT 
WHEN THE FREQUENCY RESPONSE IS GIVEN 
(2.1) Oblique Step Input 


Let » be the input variable (for example, control de- 
flection) and y be the responding variable (for example, 
pitch or rate of pitch or translational velocity of the 
aircraft). 

For the input variable n(¢) as a function of time /, 
we assume the periodic function shown in Fig. i with a 
period JT = 27/w large enough to be sure that the 
transient of the responding variable y(t) at the end of 
the half period 7/2 is practically subsided. (See 
Fig. 2.) 

We first expand the periodic function shown in Fig. 3 
in a Fourier series. Since the function n(x) is anti- 
symmetrical with respect to the point x = 0, only sine 
terms are present in the Fourier expansion: 


n(x) = . Cy sin (nx) (1) 
l 


n 
with 


n(x) sin (nx)dx (2) 


Only the odd terms of the series (1) are different from 


zero, and Eq. (2) may be written 


C, = (4 ofr sin (nx)dx, n= 1,3,... (3) 
0 
The result of the integration is 
Cn = (4b/ma)[sin (na)/n*| 
and the Fourier series (1) reads 
n(x) = (4b/ma) . [sin (na) sin (mx)]/n® (4) 
n=1,3,... 


For some purposes it is more convenient to write the 


series in complex form. Substituting 
sin (nx) = (1/21)(e"* — e~™*) 
the series (4) reads 
ee = a) 4 ° = 
n(x) = (b/ma) DO fe"? eit) fy? (5) 
n=1,3, 


~ 
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The summation is to be taken over all positive and nega 
tive odd integers n. 

In order to obtain the series for n(i), shown in Fig. 1, 
we substitute in Eq. (5) 


r T ] 
9 9 ») 


and add the mean value '/2 (compare Fig. 3 with Fig. 
Ly: 


l 1 o ernet einer — 
n(t) = = a eee p TS = ) (6) 


Tw n=1,3,... M” 


-—I1, -3 


w is the fundamental circular frequency of the periodic 
input function »(t) and of the response function y(t); 
r is the time for the change of input position (control 
deflection) from 7(0) = 0 to n(r) = 1, the rate of change 


being uniform between / = 0 and t = 7; and 7 denotes 


Vv —I. 
Eq. (6) in real form reads 


le 
a + Fe gh 


TH » , oe TN’ 


cos [nw(t — r)] — cos [nwt 


n(t) = 


( ) 
‘ 


and the summation is to be taken for positive odd in 
tegers n only. 

In order to obtain the response y(t) of the system 
to the input (/), every harmonic term in Eq. (6) has to 
be multiplied by the value of the frequency response 
function corresponding to the frequency mw. In com- 
plex notation the frequency response function, also 
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Fig. 1 (top). Periodic input (t). Fic. 2 (center). Periodic 
response y(t) to input »(t) of Fig. 1. Fic. 3 (bottom). Periodic 
function represented by Eqs. (4) and (5). 
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Fic. 5 (bott 
Periodic response y(f) to impulse input n(¢) of Fig. 4. 


Fic. 4 (top). Periodic impulse input y(t). 


called transfer function, /(inw), is the complex ampli 
tude ratio y(inw)/n(inw) = F(inw), obtained by solv- 
ing the equations of motion of the system for a steady 
state harmonic oscillation e’””. 

The transfer function is a complex rational function 
of the imaginary variable imw with real coefficients 
Therefore, F(inw) and F(—inw) are conjugate complex 
numbers, and 


ign 


for positive n | 
— ten ; 


sion \ Fe 
F(inw) = ‘ 
Fe for negative nf 


n 


F, (amplitude ratio) and ¢, (phase difference) are real 
numbers. 

Multiplying the harmonic terms of Eq. (6) with the 
corresponding values of the transfer function, 


. —inwr 
F,, ,i(nwt + gn) (< — ") 
ae 
- 


7 ia 
“ Tan=1, 3, n 
hy ¢ 


4 


Fo is the static amplitude ratio for zero frequency. 
Considering Eq. (8), the series (9) may be written in the 
real form: 


mR 9 
y(t) = > +- a Xx 
Fr § t 9 


= [nw(t — 7) +¢n] — cos (nwt + en) f 
and the summation is to be taken again for positive odd 
integers » only. The amplitude ratio F, decreases 
with 7 if nw is higher than the highest natural frequency 
of the system; the coefficient F,/n* decreases still 
faster with m; therefore, the series (9) or (10) con- 
verge rapidly once mw is beyond the highest natural fre 
quency of the system. The fundamental frequency w of 
series (9) or (10) must be chosen as a fraction of the 
lowest natural frequency and must be the smaller the 
less damped the transient is, because otherwise the solu- 
tion for a periodic input as shown in Fig. 1 would not 
be an approximation to the step input. The number 
of terms to be carried in series (9) or (10) depends on 
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EQUATIONS FOR 
the characteristics of the system. For a well-damped 
system with a not-too-wide natural frequency spectrum, 
seven to nine terms may be entirely satisfactory. The 
summation may be performed numerically or graphi- 
cally by adding vectors with the length proportional to 
F,/n°? and the direction with respect to a reference line 
given by the angles nw(f T) + g and nwt + ¢, and 
by projecting the resultant vector on the reference line. 

For 7 approaching zero (rectangular step function), 
Eq. (6) is reduced to 
© 


> (e""/in) (11 


n=1,3, ... 


n(t) = (1/2) + (1/n) 


1 


ind Eq. (9) 1s reduced to 


vs = i "eee pre , 
yl (Fo 2) rN l wr te (F,, in)e” 12 
’ 3 


2=1,3,. 
1, —3,. 


or in real form [consider Eq. (8) |: 


(F,/in) sin (nwt + 


y(t) (Fo/2) + (2/r) > 


n=1, 3,... 


[his equation was used in reference 1. 


(2.2) Comparison with Exact Solution 


If we take two discrete points on the frequency re- 
sponse curve which are close together, then the interval 
between them for » odd will be denoted 


(n + 2)iw — niw = 2iw = As 


Further, 


nw = s. Then, as w — 0, Eq. (11) becomes 


| l *—As/2 est ico e*! 
t}=s +5 | - ds +f — as| (14 
“ ee 7 =_* RY As/2 § 


We connect the two paths of integration along the 
positive and negative imaginary axes from As/2 to 
— As/2 by the right half circle around the pole of the 
integrand at s = 0. The residue of the integrand at 
’ = 0 (coefficient of 1/s in the expansion of the inte- 
grand) is one, and the value of the integral along the 


small half circle is 
fev ds = mt (15 


Combining Eqs. (14) and (15), 


+ie 
n(t) = ani) fi (e**/s) ds (16) 
~ie 
It is to be understood that the path of integration 
along the imaginary axis passes the pole of the integrand 
ats = O at the right. Eq. (16) is the well-known in- 
verse Laplace transform of the rectangular unit step 
function n(t). 
Assuming that the complex transfer function F(s 
has no poles on the imaginary axis or to the right of this 
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axis (stable systems) and assuming further that, for 
s = 0, F(O) has the real value Fp (static deflection), one 
obtains, by the same method indicated above for w ap- 
proaching zero instead of Eq. (12), the following equa- 
tion: 


* ++ §o 
yt) = 2ri) f [ F(s)/s]e*' ds (17) 


This is the inverse Laplace transform of the response 
y(t) caused by a rectangular unit step input. 

Writing again 2% = As and niw = s, Eq. (6) reads, 
for w approaching zero, 


1 1 As/2 e* | _ iad 
nf) = - + —— <i 
2 2m ia s* T 


. 


b PO ss | eS CPi es l 


results again in a residue one, and the value of the in- 
tegral along the small right half circle around the pole 
ats = Ois again 77; therefore, Eq. (18) may be written: 


oF 


1 (tiv e*1 —e 
ie / if est. EF (19) 
oat J ja, 3 T 


By the same method and assuming the above-mentioned 
properties of the complex transfer function F(s), Eq. 
(9) may be written, for w approaching zero in the form: 


1 ftie 1 —e-* 
F(s) — ————ds_ (20) 
Sad 2 
2at J —i« s? T 


Eq. (19) is the inverse Laplace transform of the oblique 
unit step function, as may be easily verified by evalu- 


rh. po 
-e "dt + | e “dl 
T 1 


e **)/rs? (21) 


y(t) = 


ating the integral 


L£[n(t)]= f n(tye * dt / 
0 J 0 


with the result 


L[n(t)] = 


Eq. (20) is the inverse Laplace transform of the re- 
sponse caused by an oblique unit step input. 

Of course, the approximate formulas [Eqs. (9) and 
(10) ] could have been derived by starting with the La- 
place transform (21) and by subdividing the integrals 
(19) and (20) in a constant part along the small half 
circle around the pole at s = 0 and in the remaining 
two parts along the positive and negative imaginary 
axes, excluding the pole, and by approximating these 
two latter integrals by an infinite series. Formulas for 
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an approximate method of calculating the transient 
from a given frequency response function, different 
from those discussed here, were derived in reference 2 
on the basis of an approximate evaluation of the inverse 
Laplace integral. Though the method of reference 2 
is more general and requires only the Laplace trans- 
form of the input function in order to be applicable to 
any given input, it seems that the use of Eqs. (9) and 
(10) or their equivalent has the advantage of a closer 
connection to a physical picture. The basis of Eqs. (9) 
and (10) is the substitution of the oblique step input 
by a periodic input, Fig. 1, and it is easy to visualize 
the conditions for good or bad approximation or con- 
vergence of the method without needing considerable 
experience with the method. 


(2.3) Impulse Input 


For the input variable n(¢) as a function of time ¢, we 
now assume the periodic function shown in Fig. 4. The 
response is Shown in Fig. 5 

The Fourier series for the input function in Fig. 4 is 
obtained by subtracting the equation 
tiw(t rT) 


in| (lla) 


n(t) = (1/2) + (1/e a le 


from Eq. (11), with the result: 


nt) = (1/7) dO (ee — &™-)/in (22 
n=1,3 
Sade 
The periodic unit impulse is obtained from this equa- 
tion by multiplication with a factor A and by letting 
the time + approach zero but keeping the product Ar 
equal to one. 
Multiplication of Eq. (22) with A = 1/7 and letting 
rt — 0 results in an indeterminate form. After apply- 
ing L’Hospital’s rule, one obtains 


@o 


n(t) = (@/r) DY el (23) 
n Ot 
1, —3, 
or in real notation 
n(t) = (2w/r) D> cos (nwt) (24) 
n=I1, 3, ... 


The response y(/) is obtained by multiplying every term 
in the series (23) with the transfer function F(inw) = 
- Feit: 


(w, TT) 2s 


y(t) = 
or, considering Eq. (8), in real notation, 


i. @ 
== Dv 


T «=i, 3, ... 


F, cos (nwt + o,) (26) 


The transition to the case where w is approaching zero 
is again made by substituting 207 = As and nwi = s. 
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Since the integrand has no pole on the imaginary axis) 
one obtains directly, instead of Eqs. (23) and (25), th 


following equations: 


~ 4:00 
n(t) = (1 2ri) | evds (27 
1 
mda 
y(t) = (I 2ni) f F(s)e*'ds (28 


Ey. (27) is the well-known inverse Laplace transform 


for a unit impulse. For a unit input during a finite 


time 7, as shown in Fig. 4, Eq. (25) reads 
- Pa sned+-on) ~én 
vit) = : ee ee [1 —e "| (25a) 
WT n=1,3, in 
1,—3 


or in real notation, considering Eq. (8), 


~ 


v(t) = — Zz. {sin (nwt + ¢n) 


wT. te . 
sin [nw(t — 7) + onl} (26a 


3) APPROXIMATE CALCULATION OF THE FREQUENCY 
RESPONSE WHEN A TRANSIENT IS GIVEN 


(3.1) Oblique Step Input 


The given response to an oblique displacement may 
be represented by the curve in Fig. 6 between ¢ = 0 and 
t = 7,2. The response curve is extended beyond the 
time / T’/2 as a periodic function with the period T. 
This period must be chosen so that at the end of the 
half period the transient is practically subsided. The 
response curve y(f) is now approximated by adding K 
periodic functions consisting of elementary trapezoids 
of the kind shown in Fig. 6. 

Each of these elementary periodic functions is ex- 
panded in a Fourier series according to Eq. (6), with 
(=f, The result is the following sum of 


AK Fourier series: 


instead of ¢. 


> l | ” eas ties 
wit) = 2) Any] s+ Zz eR 
1 


ys TW) n=1,3 n- 


29 


Eq. (9) gives the same curve y(t) in terms of the trans- 
Comparing Eqs. (9) and (29), 


” 


fer function F(inw). 
the following expressions for the transfer function are 


found: 


- 
Fe= >) Any (30) 
m=1 


—inw 


-~ K tee al 1 
+ fon U — inwt ‘ 
Fe =a Ds Aaye (— — ) 31) 
"—I1m=1 . 


For each frequency nw in the interesting frequency 


range, the summation of the A complex numbers on 
the right side of Eq. (31) must be performed. 





In} 
sti 
A, 


Th 


for 


TV axis» 


25), the 


nsf rm 
a finite 


(26a 


JUENCY 


it may 
- ( and 
nd the 
riod T, 
of the 

The 
ling K 


ez01dS 


iS @X- 
), with 
sum of 


(29 


trans- 
1 (29), 
on are 


juency 


ers on 





EQUATIONS FOR 




















ae 3 x — t 

2 w 
yt) 
~~ : fm, 
amy | 4 
tm o T | NJ 
Se ae See 
 -— a —<_ : 





| -.  —— 

i 2° w& t 
Fic. 6 (top). Response y(t) with mth element of series (29 
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Fic. 8 (bottem). Response y(t) with mth element of series (43 


If all the 7,, are approaching zero, Eq. (29) is re- 


duced to 


K x imw(t— im) 
, | | e w( m 
w(t) = sny| +- Db (32) 
1 20 wrn=13, nm 


m 


—1,—3,... 
In this case f¢,, is better chosen as the midpoint of each 
For 7,, —~ 0, Eq. (31) is re- 


time interval—see Fig. 7. 


duced to 


A 
eo NTW _— 
Fem = > Anye~™* (33 
l 


1(e asia ] )m 


For a rectangular step input, 7 > 0, Eq. (31) is reduced 
— ihwTm 


to 
. kK 
—tnwlm e l 
Dy Anye (34) 
iw m l Tm 


If the elementary steps are assumed to be rectangular, 


__ — 
re = 


Tn > 0 (see Fig. 7), Eq. (34) is further reduced to 
K 
Fev? =)> awe (35) 
m l 
For equal time intervals ¢,, 4, fm, this equation is 
identical with that used in reference 1. 
The formulas that use elementary trapezoids allow 
a smoother approximation of the given response curve 
than those that use elementary rectangles, and, there 
fore, less elements are necessary for the approximation. 
In order to obtain, in the case of rectangular unit step 
input, the connection with the exact solution, we sub- 
stitute in Eq. (35) inw = 
4,v>Oand T-— o: 


s and approach the limits 


F(s) = / (dy/dt)e-"‘dt 
0 


This equation indicates that F(s) is the Laplace trans- 
form of dy/dt. 


(35a) 


Since the system is at rest before the 


LINEAR SYSTEMS 727 
unit step input, Eq. (35a) may also be written: 


. 


F(s)/s - | vitje s¢ dt (35b) 
J0 


which is the inverse to Eq. (17) and which is the well- 
known Laplace transform of the response y(t) caused 
by a rectangular unit step input. Going to the same 
limits in Eq. (33) for the oblique unit step input, we ob- 


TS J P dy 
; e 
" Sa dt 


| = e T «x 
F(s) - / y(tje—*! dt 
Ts? 0 


his is the inverse of Eq. (20), and it was shown in Sec- 
tion (2.2) that Eq. (33b) represents the correct La- 


tain 


di (33a) 


Fis) = 
or 


(33b) 


place transform of the response y(t) caused by an ob- 
lique unit step input. 
(3.2) Impulse Input 

We assume again that the response to a periodic 
impulse input, Fig. 4, is approximated by the sum of A 
elementary trapezoids, Fig. 6, and is represented by Eq. 
(29). Comparison of Eq. (29) with Eq. (25a) leads to 
the following equation for the transfer function: 


. AK — hwrm 
a: 1 eee | - | 
Fem = ———— > Anye — 
1 Tm 


nw(1— e ‘"*" )m 
(36) 


If the periodic response curve is approximated by the 
sum of AK elementary rectangles, Fig. 7, all the 7,, in 
Eq. (36) are approaching zero and one obtains 

l 


l—e 


+ ign inwlm 27 

Fe “= inwr y% Amye (37) 
m=1 

In case of a unit impulse of infinitely small duration, 

Eqs. (36) and (37), by multiplication with 7 and letting 


7 approach zero, are reduced to 


l kK e —Mwtm __ | 
Fe?" = p Anye nwt ( ) (38) 
1 


nN ah =m Tm 
l A 


F — = > Amye~ 
l 


, inw m 


nwlm 


(39) 


In order to obtain the connection with the exact 
solution for a unit impulse, we substitute in Eq. (39) 
>Oand T 


inw = s and approach the limits A,y > @: 


(1 »f (dy/dt)e—*‘dt 
0 


This equation indicates that s/(s) is the Laplace trans- 
form of dy/dt, or, since the system is at rest before the 


F(s) = (40) 


impulse is applied, 


F(s) -f y(the—*"dt 
0 


(41) 
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known fact that the transfer function of a system is 
equal to the Laplace transform of the transient caused 
by a unit impulse. Eqs. (36) and (38), which are 
based on the approximation of the periodic response 
curve y(t) by trapezoids, allow less elements than Eqs. 
(37) and (39) because of the smoother approximation. 

Finally, a third set of approximate formulas may be 
obtained by using vertical strips instead of horizontal 
strips for the approximation of the given response curve 
y(t), see Fig. 8. 

Each of these elements is expanded in a Fourier series 
according to Eq. (23), with f t,, instead of f. The re- 


sult is: 
K - 
y(t) = b 2 Vn A,,1(w T ) =. ene —tm) 42) 
m=1 o21,8,... / 
—1,—3,... 


Comparing Eq. (42) with Eq. (25a), 
K 
ton nw inwtr ‘ 
Fy" = ———— DD ymAnte™ (43) 
l—e m=) 
This equation is valid for the case of a unit input during 
atime r. In the limit for an infinitely short time 7, but 
assuming a unit input impulse, Eq. (43) is reduced to 


Fe = > »,4,te™ (44) 
m=1 

Substituting inw = s and approaching the limits A,,t — 
0 and JT — - transforms Eq. (44) directly into Eq. (41). 
The vertical strip method may be inferior to the hori- 
zontal strip method when analyzing the response to a 
step input (Fig. 7), but it may be more convenient when 
analyzing the response to an impulse (Fig. 8). 
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* (5) APPENDIX—-SUMMARY OF FORMULAS 


(5.1) Approximate Calculation of Transients for a Given 
Frequency Response 


(5.1.1) Oblique Unit Step Input. 
Fy 2 ~ F, 
tT G=-—-+— 2 —*x 
4 < TWT n=1,3,... N° 
{cos [nw(t T) + ¢,]— cos (nwt +¢,) } 
(5.1.2) Rectangular Unit Step Input. 
a” vt) = (Fo/2) + (2/n) > (F,/n)X 
n=1,3,. 


sin (nwt + ¢,) 
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(5.1.3) Temporary Unit Input. 
2 F, 
1th wO=- LD xX 
|» Wn=1,3,...7 
=o cM {sin (nwt + gp) sin [mw(t — r) + n]} 
35.1.4) Unit Impulse Input. 


h 


1 y(t) = Qe/r) ¥ 


n=1,3,... 


F,, cos. (nwt + ¢,) 


(5.2) Approximate Calculation of the Frequency Response 

for a Given Transient 

Modes of approximation are: 

A) The given response function is approximated 

by horizontal elementary trapezoidal strips. 

(B) The given response function is approximated 
by horizontal elementary rectangular strips. 

(C) The given response function is approximated 
by vertical elementary rectangular strips. 


5.2.1.) Oblique Unit Step Input.— 
P k 
A) Fe?" = (1 'An) } A malimve — 
m=1 
k 
B) | A tt a (nw/iA p) > Anve” tNwtrm 
m=1 
5.2.2.) Rectangular Unit Step Input.— 


Ls 
\ — . — inwtm 
(A) F,e'* = (i/nw) >> AmnAnye~™ 


m=1 


k 
(B) Fe"=)> Awe" 


m=1 


5.2.3.) Temporary Unit Input.— 


k 
A) F,e*" = (r/inw An) >> AmnAmye 
m=1 


k 
Ait iw 


m=1 


(B) F,e*" = (r 


NWT “ inwtm 
> ym Ante 


~ A ee 


C) F,e* = 
Unit Impulse Input. 


. 
(A) Fye*™ = (1/n%w*) > Amn dgye™ i” 
m=1 
k 
(B) F,e'*" = (1/inw) > Anye™ 
m=1 


k 
(C) Fe" = > er 
m=1 


Abbreviations: 
oe Si 1 e Mwtm 1 
2 — = F _ Aus 
T Tm 
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Present Status of Research on 
Boundary-Laver Control 


A. E. VON DOENHOFF*® anp L. K. LOFTIN, JR." 


National Advisory Committee for Aeronautics 


SUMMARY 


A survey has been made of the present status of research 
on boundary-layer control and its possible applications in aero- 
nautics. Although the number of possible applications of 
boundary-layer control is large, only those that have received 
the most attention recently or show the most promise of pro 
ducing useful results are considered in the present paper. The 
possible applications of boundary-layer control considered are: 

1) Reduction of profile drag by the elimination of turbulent 
separation and by increasing the relative extent of laminar flow. 

(2) Increase of the maximum lift coefficient through control 
of laminar and turbulent separation 

3) The use of suction and blowing slots near the trailing edge 
of an airfoil as a means of lateral control. 

t) The use of boundary-layer control as a means of increas 
ing the efficiency of diffusers and bends. 

5) The use of boundary-layer control to influence shock 
boundary-layer interaction at high speed. 

The possible improvements in airplane characteristics result 
ing from these applications of boundary-layer control are dis 
cussed, and the general lines of future research are indicated. 


INTRODUCTION 


> PRANDTL’S FIRST PAPER On boundary lavers, 


removal of a portion of the boundary layer or the 
injection of high energy air under the boundary layer 
ave been considered as possible means for avoiding 
boundary-layer separation. More recently, removal 
of a portion of the laminar boundary layer, either 
through discrete slots or through a permeable surface, 
has been proposed as a means of increasing the stability 
and thereby the relative extent of the laminar laver at 
high Reynolds Numbers so that the skin friction drag is 
reduced. It has been suggested that the effects on the 
potential flow field of withdrawing or ejecting small 
quantities of air near the trailing edge be used as a 
means of varying the lift of an airfoil. 

A great deal of research has been conducted on vari 
ous phases of boundary-layer control, but few results of 
this vast quantity of research have found practical 
application. The purposes of the present paper are to 
state briefly the present status of information regarding 
various types of boundary-layer control, to discuss the 
possibility of improving the characteristics of aircraft 
by boundary-layer control, and to indicate the general 
lines of future research on this subject that appear to 
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Oratory. 


offer the greatest promise of producing useful results. 

The possibility of using boundary-layer control exists 
each time that the avoidance of separation or the main- 
tenance of extensive laminar flow becomes a problem. 
Because of the large number of ways in which boundary- 
layer control can conceivably be applied, no general 
conclusions regarding the applicability of boundary- 
layer control can be drawn. The scope of the present 
paper is, therefore, limited to the consideration of pro- 
posed methods of boundary-layer control as they affect 
the individual aerodynamic characteristics of an air- 
craft. The applications considered are: 

(1) Reduction of profile drag by the elimination of 
turbulent separation and bv increasing the relative ex- 
tent of laminar flow. 

(2) Increase of the maximum lift coefficient through 
control of laminar and turbulent separation. 

(3) The use of suction and blowing slots near the rear 
portion of the airfoil as a means of lateral control. 

(4) The use of boundary-layer control as a means of 
increasing the efficiency and the range of efficient oper 
ating conditions of diffusers and bends. 

(5) The use of boundarv-laver control to influence 
shock boundary-layer interaction at high speeds and, in 
particular, to eliminate boundary-layer separation fol- 
lowing the shock. 

The use of slots represerits in all cases one method of 
applying boundary-layer control, and in such cases the 
gains resulting from the use of boundary-layer control 
may depend in large measure on the design of the slots. 
For this reason, a short discussion of the status of re- 
search on the development of efficient slots is included 


at the end of the paper. 


SYMBOLS 

A = aspect ratio (b?/S) 

= wing span 
( = chord 
x = distance along chord 
t = airfoil maximum thicknes- 
Ss = wing area 
W = airplane weight 
V = free-stream velocity 
QO = volume flow per unit span 
v/V = ratio of local velocity to free-stream velocity 


W’/S = wing loading, lbs. per sq.ft 
L/D = lift-drag ratio 
= section drag coefficient 
Cr = section lift coefficient 
Cr = wing lift coefficient 
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CQ = section flow coefficient (Q/cV 
R = Reynolds Number 
ao = section angle of attack 
max. = subscript indicating maximum 
DRAG 


As mentioned in the introduction, reductions in the 
profile drag can be achieved by boundary-layer control 
through control of the turbulent layer and through ex- 
tending the length of the laminar layer. Numerous 
design calculations have indicated the comparatively 
large improvements in airplane performance to be ex 
pected from reductions in the profile drag. The im- 
provement in performance which can be expected from 
reductions in the profile drag through the maintenance 
of extensive laminar layers is shown to be even more 
marked for jet than for propeller-driven aircraft.’ The 
effectiveness of control of the turbulent boundary layer 
as a means of reducing the profile drag is considered 


first. 


CONTROL OF THE TURBULENT LAYER 


In the absence of separation, theoretically, some re- 
duction in the net drag can be obtained by sucking the 
boundary layer into the interior of the airfoil at the 
trailing edge and discharging the air thus withdrawn at 
free-stream total pressure. It seems unlikely that any 
net gain could result from such a process because of the 
necessary losses associated with the internal flow. If, 
however, fairly extensive regions of separated flow exist, 
controlling the turbulent boundary layer in such a way 
as to eliminate separation results in substantial reduc 
tion of the profile drag even when the necessary pump 
ing power is included in the drag coefficient.” * Such 
separation occurs on airfoils of moderate thickness at 
lift coefficients approaching the maximum and on ex 
tremely thick airfoil sections throughout the entire 
range of operating lift coeflicients. The reduction in 
drag observed at the higher lift coeflicients for the thin 
sections usually comes about as a by-product of at 
tempts to improve the maximum hit by boundary-layer 
control. For wings having aspect ratios less than 10 or 
12 and airfoil sections of less than 20 per cent thickness, 
such savings in drag are particularly marked only for 
lift coefficients above those for maximum L/D and, 
hence, are of relatively minor importance. 

The use of aspect ratios of the order of 15 to 25 has 
always appeared attractive from the point of view of 
lower induced drag but has not been practical because, 
for structural reasons, the root sections of such wings 
are fairly thick so that flow separation occurs at all use 
ful lift coefficients and the associated increase in profile 
drag equals or exceeds the saving in induced drag. 
Under such circumstances, the prevention of separation 
by boundary-layer control would seem to offer the possi 
bility of realizing net drag savings and increased values 
of L/D on wings of high aspect ratio. With this 
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thought in mind an experime ital investigation has re- 
cently been made at the N.A.C.A. of the characteristics 
of three NACA 6-series airfoils of 24, 32, and 40 per cent 
thickness ratio employing boundary-layer control by 
suction to prevent separation of the turbulent boundary 
layer. The airfoils were cambered to have theoretical 
design lift coefficients of 0.4. The suction was applied 
through a single slot at or slightly behind the mid-chord 
position of the airfoilsections. In order that the results 
should correspond approximately to the limited extent 
of laminar flow that is usually obtained on airplane 
wings under practical operating conditions, the tests 
were made with the leading edges of the models rough- 
Lift 
and drag data are now available for the 24 per cent thick 


ened sufficiently to cause immediate transition. 


section” and are being prepared for publication for the 
The data for the 


three thick sections, together with those for airfoils of 


52 and 40 per cent thick sections. 


12 to 21 per cent thickness employing boundary-layer 
control,'~® are sufficient to enable designers to deter 
mine the desirability of employing boundary-layer con 
trol for the purpose of improving the characteristics of 
high-aspect-ratio wings. In order to supply detailed 
design information, however, further research will be 
necessary to determine optimum slot shapes, pitching- 
moment characteristics, and the behavior of three- 
dimensional wings employing boundary-layer control. 

In order to give some indication of the possible im 
provement in wing characteristics which can be realized 
by employing boundary-layer control on wings of high 
aspect ratio, a comparison is made in Fig. | of a group 
of wings having a taper ratio of 0.4 and a varying aspect 
ratio with and without boundary-layer control. The 
wings are composed of NACA 6-series airfoils with lead- 
ing edges roughened sufficiently to cause immediate 
transition. The 
based on the structural design criterion that the ratio of 


root section thickness ratios were 


the span to root thickness should be 35 to 1, but in no 
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ASPECT RATIO, A 
(a) WINGS ALONE (b) 0.01 Cp, ADDED 
Fic. 1. Effect upon (Z/D) maz. of controlling turbulent sepa- 
ration on the thick root sections of high-aspect-ratio wings with 
leading-edge roughness. (a) Wings alone. (b) 0.01 parasite-drag 
coefficient added. 
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RESEARCH ON 


case was the root section thickness ratio made less than 
In all cases the tip had a 12 per cent thick- 
A linear spanwise distribution of thickness 


12 per cent. 
ness ratio. 
ratio was assumed. The drag values used in computing 
the values of L, D included the wake drag plus the drag 
equivalent of the suction power. 

A\ comparison of the data for the wings shown in Fig. 
| indicates that the optimum aspect ratio for maximum 
L Disinereased from approximately 11 to 20 by the use 
of boundary-layer control with an accompanying in- 
crease‘in L/D of approximately 19 per cent. The suc- 
tion power considered in the calculations is that re- 
quired to remove the necessary quantity of air from the 
surface of the wing through the slot but does not include 
any estimate of the losses that would occur in the duct 
ing and pumping equipment of an airplane. The re 
sults for the two sets of wings given in Fig. | are, how- 
ever, strictly comparable if the net efficiency of the 
pumping system from the boundary-layer control slot 
to the discharge outlet is equal to the propulsive efh- 
ciency of the main driving unit. 

Despite the fact that the L/D of the wing with bound- 
ary-layer control is higher than that of the wing with- 
out boundary-layer control, the value of the profile- 
drag coefficient at (L/D)nar, is greater for the wing 
with boundary-layer control. Consequently, the addi- 
tion of a given increment of parasite drag coefficient 
will make the comparison more favorable for the wing 
with boundary-layer control. Fig. 1b gives a compari 
son between the two sets of wings when an arbitrary 
increment of parasite drag coefficient of 0.0100 has been 
added to each. Inspection of the data of Fig. 1b indi 
cates that the optimum aspect ratio is now slightly 
higher than 20 for the wings with boundary-layer con- 
trol and that the use of boundary-layer control gives an 
increase in (L/D) mar, of approximately 50 per cent. 

It should be noted that the maximum lift coefficients 
of the thick sections are high, and, as is shown subse- 
quently, it is on wings of high aspect ratio that high 
maximum lift coefficients can be most effectively em- 
ployed. 

This application of boundary-layer control would be 
of primary interest for relatively low-speed airplanes 
where range is of the greatest importance. For ex- 
imple, the critical Mach Number obtained from low- 
speed pressure-distribution diagrams for the 40 per cent 
thick airfoil section at conditions corresponding to maxi- 
mun L/D of the wing is 0.450. In spite of a number 
of obvious difficulties, it is thought that the possibility 
of increasing this limiting speed somewhat by the use of 
a moderate amount of sweep should be investigated. 


CONTROL OF THE LAMINAR LAYER 


The first attempts to obtain reductions in the profile 
drag by increasing the relative extent of laminar flow 
consisted of the design of new airfoil shapes having the 
position of minimum pressure far back along the sur- 
face. The rearward practical limit of the position of 
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minimum pressure was dictated by the avoidance of 
turbulent separation over the rear portion of the airfoil, 
particularly for lift coeflicients outside the low-drag 
range. Three general types of boundary-layer control 
have been proposed to increase the possible relative ex 

tent of laminar flow: first, multiple slots to limit the 
growth of the boundary-layer thickness and eliminate 
laminar separation; second, special airfoil sections hay 

ing the position of minimum pressure extremely far back 
together with a single suction slot to eliminate sepa 

ration at a pressure discontinuity; third, continuous 
boundary-layer suction through a porous surface. 


The Use of Multiple Slots 


Multiple slots have two distinct and interrelated 
effects on the boundary laver; they decrease the value 
of the boundary-layer Reynolds Number and delay or 
Work on this problem has 


i] 


prevent laminar separation. 

been carried out by Holstein’ * in Germany; by Pfen- 
ninger® and Ackeret, Ras, and Pfenninger!’® in Switzer 
land; and by the N.A.C.A. in the United States. 


These investigations showed, in general, that it was pos- 


sible to extend the laminar layer in a region of adverse 
pressure gradient practically to the trailing edge with a 
small expenditure of power such that extremely large net 
drag savings were realized. The most favorable re- 
sults were obtained at a fairly low value of the Reynolds 
Number. For example, Pfenninger’s best results, 
which showed a net drag saving of 50 per cent, were ob 
tained at a Reynolds Number of 2.0 & 10°. In both 
Pfenninger’s and Holstein’s experiments, however, the 
saving in drag disappeared when an attempt was made 
to repeat the tests at higher Reynolds Numbers. 
Pfenninger attributed this adverse scale effect to in 
creasing turbulence in the tunnel as the speed was 
increased. Holstein’ was able to maintain extensive 
laminar flow up to Reynolds Numbers of 3.2 X 10°. 
However, he found that, if the slot spacing was not de 
creased as the Reynolds Number was increased, the 
power required to limit the growth of the laminar 
boundary layer would become excessive. He gave no 
explanation of his failure to achieve extensive laminar 
flow at Reynolds Numbers above 3.2 X 10°. 

An investigation is currently being carried out by the 
N.A,C.A. on a symmetrical NACA 64A010 airfoil sec 
tion of 3-ft. chord designed to maintain laminar flow to 
the trailing edge by means of suction slots up to Rey- 
nolds Numbers of the order of 25.0 & 10° The model 
is being tested at zero lift in the Langley two-dimen- 
sional low-turbulence pressure tunnel. 

The following considerations dictated the slot spac 
ing and slot size: First, the boundary-layer Reynolds 
Number should not be allowed to exceed.a definite 
value. Boundary-layer Reynolds Numbers 
sponding to transition (based on displacement thickness) 
of 6,000 to 7,500 have been measured in flight,'' and 
values of 5,000 to 6,000 have been measured on another 


wing section in the Langley two-dimensional low-turbu 


corre- 
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ic. 2. Pressure distribution and suction-slot locations for an 
NACA 64A010 airfoil section 


lence pressure tunnel.'’ In order to be reasonably con 
servative, the design value of the maximum boundary 
layer Reynolds Number for the slotted wing section was 
chosen to be approximately 2,600 at a wing Reynolds 
Number of 25.0 & 10°. Second, the spacing between 
the slots was determined from suction-power considera- 
tions. Although the suction power required to main 
tain the boundary-layer Reynolds Number decreases 
continuously with decreasing slot spacing, it was found 
that the savings of power corresponding to a slot spacing 
smaller than */, in. on the present model (a Reynolds 
Number run of 0.5 K 10°) were relatively small. Fur 
thermore, a smaller slot spacing would increase con 
siderably the mechanical difliculties of constructing the 
model. Previous investigations in the Langley low- 
turbuience tunnel indicated that the slot width should 
not be greater than the boundary-layer thickness. For 
the model under consideration this slot width was ap- 
proximately 0.005 in. 

Three parameters 
the boundary-layer Reynolds Number, the ratio of the 
slot width to the boundary-layer thickness, and the 
sufficient to 
If the values 


namely, the maximum value of 


Reynolds Number run between slots — are 
determine the design of a slot installation. 





Model of the N.A.C.A. 64A010 airfoil section equipped 
with 41 suction slots on upper and lower surfaces. 


Fic. 3. 
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of these three parameters are held constant, the slot 
spacing and slot width expressed as fractions of the 
chord will be functions of the design Reynolds Number 
For example, the slot spacing and slot width on the 
NACA 64A010 airfoil model would be 5 and 0.034 in,, 
respectively, at a Reynolds Number of 25.0 X 10° if the 
chord of the airfoil were increased to 20 ft. If, however, 
with the 20-ft. chord the design Reynolds Number were 
increased to 83 X 108, the slot spacing and slot width 
would be 1'/2 in. and 0.010 in., respectively, and the 
number of slots would increase in proportion to the 
Reynolds Number. 

The theoretical pressure distribution, together with 
the slot locations, is given in Fig. 2 for the 3-ft. chord 
NACA 64A010 airfoil model designed for a Reynolds 
Number of 25.0 XK 108. Great pains were taken in the 
construction of the model to maintain the machined 
aluminum surfaces in a smooth and fair condition. A 
photograph of the model partially disassembled is shown 
in Fig. 3. Preliminary test results indicated that not 
much difficulty was encountered in obtaining laminar 
flow over substantially the entire surface of the model 
up to a Reynolds Number of about 3.0 K 10°. As the 
Reynolds Number was increased, however, the laminar 
flow in the boundary layer became exceedingly sensitive 
to minute changes in the shape of the slot entry and 
flow quantity removed. It was found that honing the 
edges of the slot slightly with a lead pencil produced 
sufficient changes in the slot contour to affect markedly 
the maximum Reynolds Number at which laminar flow 
could be obtained over the slot. The maximum Rey- 
nolds Number at which laminar flow could be obtained 
over substantially the entire upper surface was 10.0 X 
10°, whereas the corresponding maximum Reynolds 
Number for the lower surface was 5.5 XK 10°. These 
Reynolds Numbers, although not so high as expected 
flight values, are considerably higher than those for 
which complete laminar flow was obtained in the in- 
vestigations of Holstein’ and Pfenninger® and were ob 
tained only after a great deal of effort had been ex- 
pended in trying to eliminate minute irregularities from 
the slot contours. Since the airfoil was symmetrical, 
the differences in results between the upper and lower 
surfaces are attributed to small variations in the con 
tours of individual slots. These variations were so 
slight that they could be observed only with the aid of a 
powerful magnifying glass. It was also observed that, 
once transition had occurred, no amount of suction ap 
plied downstream of the transition point restored the 
boundary layer to the laminar state. The conclusion 
drawn from this investigation is that, although the pos 
sible region of laminar flow may be extended to the 
trailing edge of an airfoil through a region of adverse 
pressure gradient at fairly high Reynolds Numbers, the 
laminar layer becomes increasingly sensitive to surface 
irregularities as the Reynolds Number is increased. 
This result is entirely consistent with those of a previous 
flight investigation'* in which no decrease in sensitivity 
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of the laminar boundary layer to surface irregularities 
was observed to result from the installation of a number 
of suction slots on a wing panel. In view of the ob- 
served increasing sensitivity of the laminar layer to sur- 
face irregularities with increasing Reynolds Numbers 
and the difficulties that have been experienced in the 
past in obtaining the design extent of laminar flow on 
low-drag airfoils on operational airplanes, the use of 
suction slots to increase the possible extent of laminar 
flow does not appear to be attractive. The practica- 
bility seems especially limited when consideration is 
given to the extreme difficulty of manufacturing and 
maintaining sufficiently accurate slot contours. 


Airfoils Designed Especially for Boundary-Layer Control 


[he second method of overcoming the limitations on 
the design extent of laminar flow imposed by the con- 
sideration of turbulent separation at the rear of the air- 
foil was suggested by Griffith and discussed in some de- 
tail by Goldstein in his Wright Brothers lecture.'4 The 
original basic idea of this method of approach was to de- 
sign an airfoil that had favorable pressure gradients 
over the entire region from leading edge to trailing edge. 
In order to obtain a closed shape consistent with this 
condition, it was necessary that the pressure increase 
discontinuously at some point along the airfoil surface, 
Suction was to be introduced at this singular point in 
order to enable the flow to follow the contour without 
separation. A typical pressure distribution and corre- 
sponding airfoil profile, taken from Goldstein’s paper, '* 
ire shown in Fig. 4. Because of the necessarily concave 
nature of the surface downstream of the pressure dis- 
continuity and the corresponding Goertler type of insta 
bility, it was not possible to obtain laminar flow down- 
stream of the suction slot except at extremely low Rey- 
nolds Numbers. Consequently, in spite of the favor- 
able pressure gradient over the rear portion of the air- 
foil, laminar flow could be expected only in the region 
upstream of the slot. Later, airfoils of this type were, 
therefore, designed with the pressure discontinuity and 
associated suction slot at a more rearward position than 
shown in Fig. 4. More rearward positions of minimum 
pressure and correspondingly lower drag coefficients 
would be. feasible with this type of airfoil section than, 
for example, with NACA 6-series sections without bound- 
iry-layer control, provided laminar flow were obtained 
up to the slot. If, however, laminar flow were not ob- 
tained up to the slot, it seems unlikely that the suction 
uirfoil would show an appreciably lower drag coefficient 
than that of a plain airfoil section designed to have 
minimum pressure at the assumed forward position of 
transition. Practical airfoils can be designed with the 
position of minimum pressure as far back as 60 per cent 
of the chord. Experience with operational airplanes 
having low-drag wings, however, indicates that laminar 
flow usually extends over a distance of no more than 15 
to 20 per cent of the chord back from the leading edge.'® 
The difficulty appears to be not only the presence of in- 
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Fic. 4. Profile shape and pressure distribution of a sym- 
metrical suction airfoil section of approximately 16 per cent 
thickness (reference 14). 


accuracies in construction but also the accumulation of 
insects and dirt associated with the necessarily exposed 
nature of the wing surfaces. There is no reason to ex- 
pect that the laminar boundary layer over the forward 
portion of suction airfoils of the Griffith type would be 
noticeably less sensitive to small surface irregularities 
than the corresponding region for NACA 6-series air 
foils without suction. Unless the certainty of obtaining 
extensive laminar flow over more than the first 60 per 
cent of the airfoil chord can be made considerably 
greater than it is at present, it is not likely that an air- 
plane designer would feel inclined to compromise the de- 
sign of his airplane to the extent of using this type of suc- 
tion airfoil. Since the advantages of extensive laminar 
flow are well known and the drag corresponding to vari- 
ous extents of laminar flow can be calculated theoreti- 
cally, further research on the design of the Griffith-type 
airfoils and on their experimental characteristics under 
ideal conditions is much less urgent than is research on 
methods of insuring the realization of extensive laminar 


flow. 


Area Suction 

A basic difficulty of obtaining laminar flow on air- 
planes is the sensitivity of the laminar boundary layer 
at high Reynolds Numbers to surface defects that are 
sufficiently small so as to be almost unavoidable. The 
only method of boundary-layer control that offers even 
any theoretical hope of reducing the sensitivity of the 
laminar layer to such small disturbances consists of con- 
tinuous suction through a porous surface. The theory 
of the stability of laminar boundary layers to small two- 
dimensional disturbances was developed by _ Toll- 
mien!® !7 and Schlichting'* and checked experimentally 
for the Blasius velocity distribution by Schubauer and 
Skramstad.'® The theory was extended by Schlichting 
and others?°-** to include the effects of variations in 
pressure gradients and the effects of blowing or sucking 
through a porous surface on the stability of the laminar 


laver. One of the most important conclusions of this 
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Model of the N.A.C.A. 64A010 airfoil section having the 
center section covered with sintered bronze 


Fic. 5 


theoretical work (dealing only with small two-dimen 
sional disturbances) is that continuous suction through 
a porous surface markedly stabilizes the laminar layer 
with respect to small disturbances and that the quan 
tity of air that has to be removed to achieve this marked 
stabilizing effect 1s extremely small. For example, the 
theory indicates that the lower critical boundary layer 
Reynolds Number based on the displacement thickness 
for the flow over a flat plate with zero pressure gradient 
is increased from a value of approximately 420 without 
suction®® to an asymptotic value of 55,000°% to 70,000, °° 
with an amount of suction corresponding to a compo 
nent of velocity normal to the plate of the order of 0.01 
of 1.0 per cent of the free-stream velocity. 

An experimental investigation to determine the effects 
of continuous suction on the drag of an NACA 64A010 
airfoil of 3-ft. chord its now being carried out by the 
N.A.C.A. in the Langley two-dimensional low-turbu 
lence tunnels. thick 
sheet of porous bronze which is made up of powder con 


The skin of the model 1s a *, y-1n. 


sisting of approximately spherical particles of such size 
that all the particles will pass through a 200-mesh screen 
but none through a 400-mesh 
wrapped in one continuous piece from the trailing edge 


screen. The sheet is 
on the upper surface around the leading edge to the 
trailing edge on the lower surface. The porous region 
has a span of | ft. situated in the center of the 3-ft. span 
model. A photograph of the model is shown in Fig. 5. 
Although some waviness was present in the model, the 
chordwise waves were actually much less severe than 
seems to be indicated in this photograph. 

Cypical results in the form of drag coefficient against 
flow coefficient are shown in Fig. 6 for a Reynolds Num 
ber of 6.0 & 10°. 
the trailing edge indicated that laminar flow was ob 


Boundary-layer surveys taken near 
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tained over virtually the entire porous surface of the 
model for flow rates corresponding to the lowest drags 
obtained. These data show that substantial net say 
ings in drag can be obtained and that completely lami 
nar flow can be maintained even when the model is not 
quite aerodynamically smooth and fair. The fact that 
the model was not quite aerodynamically smooth and 
fair is shown by the comparison of the drag coefficient 
for the boundary-layer control model with sealed sur 
face and the corresponding drag coefficient of the solid, 
fair, and smooth model of the same airfoil section. In 
sealing the porous model, the surface texture was not 
altered. At Reynolds Numbers substantially higher 
than 6.0 X 10°, it was not possible to obtain any net re 
duction of drag. This adverse scale effect appears to be 
associated with the particular sample of material used 
in the investigation. The pressure drop across the 
porous material is directly proportional to the flow veloc 
ity through it, so that the chordwise distribution of in 
flow velocity becomes increasingly nonuniform not only 
with decreasing flow coefficient but also with increasing 
Reynolds Number. The flow coefficient corresponding 
tO C4,j,, ata Reynolds Number of 6.0 X 10° is somewhat 
greater than that indicated as theoretically necessary 
with a uniform inflow velocity to obtain the desired sta 
bility. 
inflow velocity necessary to avoid local regions of out 


Relative to free-stream velocity, the minimum 


flow increased with increasing Reynolds Number. At 
high values of the flow coefficient, 1t was rather difficult 
to judge whether the boundary layer was laminar or tur 
bulent. In general, however, the results seem to indi 

cate that complete laminar flow was obtained provided 
there were no local regions of outflow over the surfac« 

At Reynolds Number much above 6.0 X 10°, the flow 
coeflicient necessary to satisfy this condition was so 
large that no net saving in drag was obtained in spite ot 
the fact that the internal structure of the model was di 
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vided into a number of compartments with separately 
adjustable suction pressures. It is planned to continue 
the investigation using a much more dense porous ma- 
terial. It should be pointed out that the pressure drop 
through the porous surface itself is sufficiently small 
compared with the free-stream dynamic pressure for the 
flow rates of interest that the pressure drop can be in- 
creased several fold without materially affecting the suc- 
tion power requirements. 

In addition to the stabilizing effect indicated in the 
discussion of the data of Fig. 6, a further indication of 
the stabilizing action was obtained by spanwise drag 
surveys in the neighborhood of the juncture between the 
porous and solid portions of the surface of the model. 
This juncture was not completely smooth. As a result, 
transition spread inward from the juncture over the 
porous region and decreased the spanwise extent of the 
low-drag region behind the model. It was noticed that 
the spanwise extent of the low-drag region increased 
with increasing inflow velocity, which indicates that 
continuous suction decreases the angle of spread of tur- 
bulence. 

Before any recommendations can be made regarding 
the use of continuous suction on airplane wings, not only 
must the feasibility of obtaining substantial reductions 
in drag be determined at higher Reynolds Number but, 
more importantly, the effects of surface irregularities 
such as are likely to be present under practical operating 
conditions must be found. 


Maximum LIF1 


Usable Maximum Lift Coefficient 


One of the earliest applications of boundary-layer 
control to receive attention was that of increasing the 
maximum lift coefficient. The gains in performance 
associated with such an improvement in wing character 
istics were thought to be obvious. It is not at all cer 
tain, however, that such is the case. For example, a 
recent analytical investigation of a conventional low 
speed airplane having a pay load of 5,000 Ibs.*" has indi- 
cated that the gains in take-off performance resulting 
from increasing the available maximum lift coefficient 
from values of the order of 5.0, which can be obtained 
without boundary-layer control, up to a value of ap- 
proximately 5.0, which can be obtained only with bound- 
ary-laver control, did not result in a proportionate de 
crease in the total take-off distance. The improvement 
in take-off performance appeared to be relatively un- 
important for aspect ratios much less than 15. The re 
sults of the analysis are consistent with results of Ger 
man flight tests of two airplanes incorporating bound- 
ary-layer control to increase maximum lift coeffi 
cient.** It should be pointed out that the take-off dis 
tance considered to be of primary importance in these 
investigations was the distance required to clear a 50-ft. 
obstacle. In nearly all cases, increases in the maximum 
lift coefficient resulted in decreases in the ground run, 
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which might be of considerable importance in special 
problems such as those encountered in the design of air- 
craft for carrier operation. 

These investigations served to clarify considerably 
current concepts regarding the usefulness of high maxi- 
mut lift coefficients for the particular take-off problem 
studied. Similar studies of both take-off and landing 
performance are badly needed for other types of air- 
craft, particularly those designed primarily for high 
speed performance and having extremely thin wings or 
wings of unusual plan form. Although the usable maxi- 
mum lift coefficients for most of the proposed high- 
speed wing configurations are probably lower than those 
of wings of more conventional plan form because of the 
associated high induced drags and low take-off thrusts, 
there does seem to be a possibility of improving the 
landing and take-off characteristics of such high-speed 
configurations by increasing the maximum lift coeffi- 
cients above the present extremely low values. There 
appear then to be two possible fields of application for 
boundary-layer control to increase the maximum lift 
coefficient: first, to relatively low-speed airplanes hav 
ing wings of extremely high aspect ratio; and, second, 
to high-speed airplanes with wings that have extremely 
low maximum lift coefficients. 


Low-Speed Configurations 


For conventional wings of high aspect ratio, methods 
exist for predicting the wing characteristics from airfoil 
section data. The discussion of methods of improving 
the maximum lift of conventional wings is, therefore, 
concerned with results that have been obtained from 
two-dimensional investigations of airfoils with bound- 
ary-layer control and other high-lift devices. 

For smooth airfoils at all reasonably high angles of 
attack, laminar separation occurs near the leading edge, 
but below the maximum lift coefficient the flow re 
attaches itself to the surface forming a turbulent bound- 
ary layer. The amount of pressure recovery that can 
occur before the turbulent boundary layer separates de- 
pends markedly on the details of the flow conditions 
associated with the initial forming of the turbulent 
boundary layer. Turbulent separation near the trail- 
ing edge and the laminar separation near the leading 
edge have a regenerative effect upon each other.” 
Maximum lift finally occurs as a result of either a pro- 
gressive forward movement of separation from the trail- 
ing edge or permanent separation of the laminar bound- 
ary layer near the leading edge. Because of the re- 
generative effect, increases of maximum lift coefficient 
on almost any given airfoil can be obtained by delaying 
either form of separation. The larger effect, however, is 
generally obtained by delaying the type of separation 
that finally results in complete flow breakdown. For 
example, the thicker airfoils with blunter leading edges, 
which have round-top lift curves, generally can be im- 
proved most by delaying separation of the turbulent 
boundary layer; whereas the largest increases in maxi- 
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inum lift of the thinner sections can be obtained by con- 
trolling separation near the leading edge. In any case, 
if boundary-layer control is used to prevent one type of 
separation, maximum lift will then be limited by the 
other type. 

The type of boundary-layer control that has received 
the most attention has been that which delays turbulent 
separation over the rear portions of airfoils of 12 per 
cent thickness and greater. Turbulent separation can 
be delayed either by removing a portion of the low- 
energy air in the boundary layer or by injecting high- 
energy air under the boundary layer. Boundary-layer 
control is effective in increasing the maximum lift co 
efficient either with or without other high-lift devices. 
Their use in connection with boundary-layer control, 
however, generally has two advantages: First, the 
values of the maximum lift that can be obtained are 
greatly increased; and, second, the angles of attack for 
maximum lift are not excessive when trailing-edge high- 
lift devices are employed. A comparison of the most 
common mnethods of controlling the turbulent boundary 
layer is given in Fig. 7. The figure shows a plot of 
maximum lift coefficient as a function of blower power 
for a given wing loading. The data presented were ob 
tained in the United States,’ Germany,* and England.?*! 
The choice of the most effective method of boundary 
layer control is seen to depend upon the power expendi 
ture per unit wing area. The data indicate that for the 
lowest power expenditures the mid-chord suction slot, in 
combination with a trailing-edge double-slotted flap, is 
most effective. Extremely high maximum lift coeffi 
cients can be obtained with an arrangement whereby 
air is blown over the flap— but only with relatively large 
expenditure of power. The arrangement whereby air 
is withdrawn in the neighborhood of the flap hinge may 
be slightly better than the other two arrangements for 
intermediate power expenditures. 

Some of the results? ‘~® of a systematic investigation 
of boundary-layer control on smooth airfoils of various 


thickness ratios are given in Fig.8. In each case bound- 
ary-layer control was applied through a single suction 
slot located at the approximate mid-chord position 
The increment of maximum lift coefficient due to bound- 
ary-layer control increased progressively with airfoil 
thickness ratio. The reason for the relatively small 
increments in maximum lift observed for the thinner 
sections is that for these airfoils maximum lift was origi 
nally limited by permanent laminar separation near the 
leading edge. In all cases with suction applied, max- 
imum lift finally occurred as a result of permanent 
laminar separation near the leading edge. 

An obvious method of further increasing the maxi- 
mum lift coefficient is to delay or eliminate leading-edge 
separation. This can be done by the use of leading- 
edge slats or flaps or by the use of boundary-layer con- 
trol. The effect of the addition of a leading-edge slat 
to the 12 per cent thick airfoil with boundary-layer con 
trol and double-slotted flap‘ can be seen in Fig.9. Sub 
stantial increments in maximum lift are gained by the 
use of the leading-edge slat such that maximum lift co- 
efficients of the order of 4.0 are possible for all of the air- 
foils of 12 to 24 per cent thickness. 

In the hope that some form of boundary-layer control 
might be more effective or convenient in controlling 
leading-edge separation than slats or flaps, several in 
vestigations have been made. The types of boundary 
layer control investigated include the location of slots 
near the leading edge and the use of a porous leading 
edge.'+ 32, 83 All data that are available from these in 
vestigations are for airfoils employing leading-edge 
boundary-layer control alone without means for con 
trolling separation over the rear of the airfoil. As might 
have been expected, the boundary-layer control elimi 
nated leading-edge separation, but turbulent separation 
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boundary-layer contro! by suction and double-slotted fla; 
R = 2.2 X 10 





bound- 
suction 
sition, 
20und- 
airfoil 
small 
hinner 
} OFigi- 
ar the 
, Max- 
lanent 


maxt- 
r-ecl ge 
iding- 
r con- 
e slat 
rcon 
Sub 
y the 
ft co- 
e air- 


ntrol 
ling 
ul in 
lary 
slots 
ding 
e in- 
edge 
con 
ight 
imi 
tion 


lol 





RESEARCH ON BOUND 





BLC- Cg = 015 


4.05 
lll Lim ULE 


3.75 


3.1 Limzzzee Le 





R = 6.0x10® 








Oo 

Fic. 9. Influence of a leading-edge slat upon the maximum 

lift obtained for a 12 per cent thick airfoil equipped with a 

double-slotted flap and boundary-layer control by suction at the 
mid-chord position. 


over the rear of the airfoil limited the maximuin lift to 
values of the order of those obtainable with a slat. 
Further research is needed in order to determine whether 
boundary-layer control applied to the leading edge 
of the thinner sections will prove more effective than 
leading-edge slats when used in conjunction with other 
types of boundary-layer control and high-lift devices. 
Boundary-layer control by continuous suction near the 
leading edge may have some advantages over discrete 
slots or leading-edge slats in that, presumably, detailed 
investigation of individual sections would not be neces- 
sary to obtain optimum configurations. 


High-Speed Configurations 


Wing configurations that have been designed pri 
marily to obtain good aerodynamic characteristics at 
high Mach Numbers generally have airfoil section 
thickness ratios of less than 12 per cent and may have 
considerable amounts of sweep. Both of these char 
acteristics lead to low values of the maximum lift co 
efficient. The low maximum lift of the thin sections is 
caused by relatively early separation of the flow from 
the leading edge. The largest improvements in the 
maximum lift would therefore’ be expected to occur as a 
result of control of the flow separation near the leading 
edge. Recent investigations at the N.A.C.A. Langley 
Aeronautical Laboratory have shown that the use of a 
plain, drooped leading-edge flap in conjunction with a 
plain trailing-edge flap increased the maximum lift of a 
6 per cent thick airfoil section from 0.78 to 1.89. At 
least equally large increments in the section maximum 
lift coefficient could probably be obtained by substitut- 
ing boundary-layer control for the flap at the leading 
edge, but the pressure difference through which the 
boundary-layer control blower would have to operate 
would be extremely large. This pressure difference 
would probably be a substantial fraction of the absolute 
pressure with normal landing speeds for airfoils of the 
order of 6 per cent thickness. It is questionable whether 
this application of houndary-layer control would be 
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sufficiently more effective than the simple leading-edge 
flap to warrant its use. No final conclusion can be 
reached, however, until data are obtained on the pres- 
sure and flow-quantity requirements. 

The maximum lift coefficients of the swept-type wings 
now being used for high-speed aircraft are generally ex- 
tremely low. Although at the present time there*are 
not generally available any data on the application of 
boundary-layer control to swept wings, the flow phenom- 
ena believed to result in the occurrence of maximum 
lift on swept wings are briefly discussed in order to indi- 
cate by what means the maximum lift of such wings 
might be improved by the use of boundary-layer con- 
trol. According to the concepts of simple sweep theory, 
the characteristics of individual sections of an infinite 
yawed wing depend upon the component of velocity nor- 
mal to the leading edge. The characteristics of finite 
sweptback wings are, however, rather strongly in- 
fluenced by three-dimensional effects not present in the 
case of the infinite yawed wing. 

The distribution of shed vorticity has two important 
adverse effects upon the characteristics of the swept- 
back wing: First, the induced vertical velocity field 
shifts the spanwise center of pressure outboard as the 
sweepback is increased; and, second, an effective nega 
tive camber is induced in the sections near the tip. The 
resultant effect upon the flow is that pressure peaks 
near the leading edge of the outboard sections tend to be 
accentuated. Asa result, the tip sections of sweptback 
wings would be expected to stall sooner than do those 
near the root, with the stall originating near the leading 
edge. Early tip stalling is further provoked by the fact 
that the spanwise pressure gradient exercises a measure 
of boundary-layer control on the sections near the root 
and by the fact that the distribution of shed vorticity 
induces an effective positive camber in these root sec 
tions. Consequently, the first steps in attempting to 
improve the low-speed characteristics of such wings 
should be the delay of leading-edge separation on the 


outboard portions of the wing. 


LATERAL CONTROL 


The effectiveness of boundary-layer control as a lat 
eral-control device depends upon the sensitivity of the 
lift of an airfoil section to details of the flow conditions 
at the trailing edge. Various investigations have been 
made in Germany“ of the effect upon the direction of 
the stream line leaving the trailing edge—and, conse- 
quently, the lift-—of discharging and withdrawing air 
through different arrangements of slots located at or 
near the trailing edge. Several of these devices proved 
quite effective in changing the lift coefficient at a given 
angle of attack in much the same way as an aileron acts. 
Although there is little reason to doubt the effective- 
ness of these devices, at least at subcritical speeds, or 
the fact that they might lead to extremely light control 
forces, it is not evident that they would prove to be 
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simpler or more reliable. than conventional lateral- 
control devices with boosters where necessary. 


DIFFUSERS AND BENDS 


Efficient diffusers are even more important on jet 
type airplanes than on airplanes with conventional 
power plants, because any losses in the diffuser would 
not only represent an increment of drag but would also 
greatly decrease the efficiency and output of the jet 
engine itself. It is extremely difficult to determine any 
general rules or design criterions for the use of boundary- 
layer control on diffusers because of the marked effect 
of the initial conditions of the boundary layer at the 
entrance of the diffuser on the pressure-recovery char- 
acteristics and the rapidity with which the diffusion can 
be accomplished without encountering serious losses. 
In many cases, the entrance to the duct is situated well 
back on a body—as, for example, a scoop inlet on a 
fuselage. In these cases, improvements in the efficiency 
of the diffuser can be obtained by removing a portion of 
the boundary-layer air ahead of the duct as well as 
within the duct itself. The bypassing of low-energy 
boundary-layer air from the duct entrance is now gen- 
eral practice for nearly all scoop-type air inlets. The 
use of suction to avoid separation within a diffuser is 
not a new principal. Unfortunately, however, general- 
ized data giving the quantity and pressure require 
ments are not available. Such data for several spe- 
cific configurations, however, are given in an N.A.C.A. 
paper by Biebel.* In general, it appears that high 
efficiency diffusion can be obtained with the use of a 
single suction slot withdrawing a quantity of air of the 
order of 5 per cent of the total quantity of air passing 
through the duct. In many cases, the improvement in 
airplane performance gained as a result of the improved 
efficiency of diffusion might more than counterbalance 
the losses associated with withdrawing the required 
boundary-layer air. Because, in general, the static 
pressure within the diffuser is higher than free-stream 
static pressure, no auxiliary pumping equipment is 
necessary. This application of boundary-layer control 
appears attractive and should be considered whenever 
the problem of efficient diffusion in a short distance 
arises. Because of the varied nature of individual ap- 
plications, however, it is difficult to outline a systematic 
research program that would provide adequate data. 
Future research would probably most profitably deal 
with proposed specific installations. — 

The application of suction to prevent separation in 
bends does not appear so attractive as that just dis- 
cussed, because the local pressure at the point where 
boundary-layer control is required is generally fairly 
low as compared with free-stream static pressure. 
Furthermore, considerable improvement in the effi- 
ciency of bends can be obtained by the proper use of 
guide vanes. 

Little detailed information is available on the use of 
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blowing slots to improve the flow in diffusers and bends, 
One such installation has been made, however, in the 
exit cone of the Langley high-speed 7- by 10-ft. wind 
tunnel where a comparatively small amount of air, hay- 
ing a total pressure equal to that in the center of the 
tunnel, is introduced under the boundary layer in the 
exit cone. It was found that these blowing slots had 
little effect on the energy ratio of the tunnel, but they 
eliminated the unsteadiness of the flow in the tunnel, 
It is possible that a corresponding arrangement in the 
entrance diffuser of a jet engine installation might have 
a beneficial influence on the steadiness of the flow enter- 


ing the compressor. 


BOUNDARY-LAYER CONTROL AT HIGH MACH NUMBERS 


Another possible application of boundary-layer con- 
trol is the control of separation following a shock at high 
Mach Numbers. It is fairly well agreed that one of the 
principal reasons for the rapid increase in drag above 
the critical Mach Number is the separation of the flow 
from the surface that accompanies shock formation 
rather than the losses in the shock itself, at least at low 
supercritical Mach Numbers. It also seems likely that 
the position of the shock is strongly affected by bound- 
ary-laver conditions. Boundary-layer control as a 
method for preventing separation following a shock has 
been investigated to some extent in Germany*’ and 
England.*’ These tests indicate that the external drag 
can be reduced at some Mach Numbers, but in many 
cases the power requirements were about equivalent to 
the saving in wake drag. These investigations, how 
ever, are not necessarily conclusive. In each case bound- 
ary-laver control was applied through a single suction 
slot. If the primary purpose of the boundary-layer 
control is to eliminate the flow separation associated 
with shock formation, there is some doubt as to the 
effectiveness of any single slot configuration because of 
the large variation of the shock position with Mach 
Number and artgle of attack. 
in drag coefficient and elimination of buffeting and the 


The possible reductions 


fact that boundary-layer control may tend to reduce or 
at least postpone to higher Mach Numbers the large 
erratic changes in lift and pitching moment makes fur- 
ther investigations of boundary-layer control at high 
Mach Numbers seem extremely important. These in- 
vestigations should include not only the effects of indi- 
vidual suction and blowing slots but also the effects of 
suction through a porous surface. The purpose of the 
porous surface in this case would not be so much to 
maintain extensive laminar flow as to ensure that suc 
tion would always be applied in the vicinity of the shock. 

The application of boundary-layer control at super- 
sonic speeds is not clear. On the one hand, analysis 
indicates that for bodies of optimum fineness ratio the 
skin friction accounts for one-half to two-thirds of the 
total drag. The investigations of Theodorsen and 
Regier** indicate that the laws for turbulent-boundary- 
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laver skin friction are not greatly different at supersonic 
speeds than at subsonic speeds. Reductions of the 
skin friction drag must, therefore, as at subsonic speeds, 
come about through an increase in the relative extent of 
laminar flow. On the other hand, the details of the 
shock formation at the trailing edge of a supersonic air- 
foil section appear to have a large effect on the drag. 
In general, thickening of the boundary layer near the 
trailing edge appears to increase the trailing-edge pres- 
sure and thereby results in a decrease of the pressure 
drag.** Boundary-layer control or any other effect 
that would tend to increase the extent of laminar flow 
and decrease the skin friction drag would appear, there 
fore, to have a tendency to increase the pressure drag. 
Hence, boundary-layer control at supersonic speeds pre- 
sents a much more complex problem than at subsonic 
speeds and is badly in need of thorough investigation. 


SUCTION-SLOT DESIGN 


Che requirements for a good suction slot depend upon 
whether the slot is to be used for removing a portion of 
the laminar or turbulent boundary layer. The primary 
requirement for slots designed to extend laminar flow is 
that the slots themselves shall not introduce any dis 
turbance that will cause transition. Slots have been 
designed that satisfy this condition in several investi- 
gations,’~* '* and further work on the design of suc 
tion slots for laminar layers is discussed by Fage and 
Sargent.” Losses in slots designed to control the lami 
nar boundary layer are usually not of critical impor- 
tance because the amount of air withdrawn at a single 
slot is small in a correctly designed installation. Fur- 
thermore, the velocity with which the air is withdrawn 
cannot be very large without disturbing the laminar 
layer and causing transition. It is necessary, of course, 
that the flow into the slot be stable, and this question 
has been investigated by Smith and Roberts. ' 

The conditions affecting the design of slots to operate 
ina turbulent boundary layer are quite different. fn 
this case the external flow is relatively insensitive to de- 
tailed changes in the slot design and is affected pri 
marily by the quantity of air withdrawn. Changes in 
the design of the suction slot do, however, have a 
marked effect on the internal losses. 
more important than for laminar layers because of the 


These losses are 


relatively large quantities of air withdrawn in each slot. 
Investigations to develop efficient slot configurations for 
turbulent boundary layers have been carried out by 
Pierpont'! and Gerber’. The problem of reducing the 
losses following entry of the air into the slot is primarily 
that of designing efficient diffusers and bends. This 
problem has been considered in some detail by Raw- 


cliffe, 4? 


CONCLUSIONS 
*. 
Consideration of the present status of research on 
various methods of boundary-layer control indicates the 


following conclusions: 
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(1) For relatively low-speed, long-range aircraft, 
boundary-layer control may be effectively employed to 
eliminate turbulent separation on thick root-section air- 
foils such that wings of higher aspect ratio may be em 
ployed to give improved values of the lift-to-drag ratio. 
The data on which the analysis was based were obtained 
for airfoils having completely turbulent boundary 
layers extending back from the leading edge. 

(2) In order to obtain extensive regions of laminar 
flow and correspondingly low profile-drag coefficients 
such as might be obtainable with NACA 6-series airfoils 
or airfoils of the Griffith type, some means must be 
found fordecreasing the sensitivity of the laminar bound 
ary layer to surface imperfections that are apt to occur 
under practical operating conditions. The use of mul 
tiple slots does not appear to decrease the sensitivity of 
the laminar boundary layer. Although information 
regarding the effects of area suction is not sufficiently 
complete to be conclusive, such data as are available 
indicate that area suction does have some stabilizing 
action and that the suction power requirements are 
small. Further research should be carried out on 
boundary-layer control by area suction. 

(3) By the appropriate use of boundary-layer control 
maximum lift coefficients of the order of 4.0 can be ob 
tained for airfoil sections of 12 per cent thickness and 
above without the expenditure of excessive amounts of 
power. Maximum lift coefficients of the order of 3.0 
or 4.0 are effective in decreasing the take-off distance 
only for airplanes having wings of extremely high aspect 
ratio. 

(4) The use of boundary-layer control as a means of 
increasing or decreasing the lift of an airfoil independ 
ently of the angle of attack has been the subject of sev 
eral investigations. These methods proved to be effec 
tive. It is not clear, however, that they would be any 
simpler or more effective than conventional control sur- 
faces. 

(5) The use of suction to control separation appears 
to be a particularly convenient method of increasing the 
efficiency of short diffusers because the pressure at the 
suction slot usually is sufficiently high to eliminate the 
need for auxiliary pumping equipment. Because of the 
varied nature of individual applications, however, future 
research on this problem would probably most profit- 
ably deal with proposed specific installations. 

(6) Several short investigations have been made of 
the effect of boundary-layer control on the drag at 
supercritical Mach Numbers. These tests showed little 
net decrease in the drag. The method employed for 
applying the boundary-layer control, however, did not 
appear to be the optimum. Further research on bound 
ary-layer control at supercritical speeds is necessary in 
order to explore more completely the possibilities of the 
application of boundary-layer control not only for re 
ducing the drag but also for improving the lift and mo 


ment characteristics. In particular, it is felt that the 
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use of continuous or area suction should be investigated 
at high Mach Numbers. 
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Heat Transter in Laminar Compressible 
Boundary Layer on a Porous Flat 
Plate with Fluid Injection 


SHAO WEN YUAN{T 
Polytechnic Institute of Brooklyn 


SUMMAR\ 


A theoretical investigation was made of flow of hot fluid over 
, porous flat plate under the condition of uniform fluid injection 
The von Karm4n integral rela 
tion of the momentum equation was derived with the uniform 
velocity of injection taken into consideration. The velocity pro- 
file was assumed as a polynomial of the fourth degree. 

Solution of the above equation gave the relationship of length 
in the direction of flow to the boundary-layer thickness. The 
velocity profiles and the temperature profiles were then calculated 
fora Prandtl Number equal to unity 

rhe relation between wall temperature and the rate of coolant 
injection were calculated for different Mach Numbers, Reynolds 
Numbers, and viscosity variation with temperature. 


NOMENCLATURI 


\ = distance along the plate in the direction of the 
free stream measured from stagnation point 

y = distance perpendicular to the plate 

7 = velocity at any point along the plate 


= velocity at any point perpendicular to the plate 
= constant injection velocity fluid 


= velocity outside the boundary layer 
= thickness of boundary layer 
\ = (195/6%_) — T; = nondimensional boundary-layer 
thickness 


= (Ux/) (povo/p: l nondimensional distance 
along the plate 


= specific heat at constant pressure 


t = coefficient of heat conduction 
= dynamic viscosity at wall and outer boundary 
layer 
Yo Piy Po = mass fluid density at wall and outer boundary 


layer, and of coolant 
ia = kinematic viscosity of fluid at wall and outer 
boundary layer 
= fluid temperature 
= fluid temperature at outer boundary laver 


= temperature of the coolant to be injected 
(71 — Tx)/T 
) se PWle 
= y/é 
R = UL/», Reynolds Number 
L = length of the plate 


1 
1 
Te = temperature of the plate 
I 
1 
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his paper is based on work supported by Project Squid, U.S 
Navy. The author’s thanks are due to his colleagues, C. Chin 
and N. Ness, for checking and carrying out the derivation and 
numerical work. 
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M = speed ratio or Mach Number of the free stream 
= specific heat ratio 
' ur — 1 mpk—1._. 
r, ogee ee mae M? for p ~ WT 
nr — 1 npk — 1 
Ty ean + ~ —— M?for uw ~ T°*’4 
8 8 2 
oe = (my — 1)/12 forp ~ VT 
= (np — 1)/8forn ~ T*’* 
Q = povo 
w = pi U 


INTRODUCTION 


| peor OF COMBUSTION CHAMBERS in rocket and jet 
motors, which are subject to extremely high tem- 
perature, by means of injecting coolant through a po 
rous wall is a new and effective cooling method. Thisis 
the so-called sweat-cooling method. The combustion 
chamber walls that are to be cooled are built from po- 
rous materials, and a cold medium, called coolant in the 
form of a gas or liquid, is pressed through the pores. 
On emerging from the wall, the coolant first picks up 
heat from the wall and lowers the wall temperature, 
and then it forms a protective layer on the surface 
exposed to the gas, which diminishes the heat transfer 
from the hot gas to the wall. Since the magnitude of 
the boundary-layer thickness along the combustion 
chamber liner is known to be comparatively small, the 
problem can be treated by investigating the boundary 
laver along a flat plate under the same flow condition. 
It is believed that this approximation will yield some 
indication of the actual case. 

The theoretical study of the aerothermodynamic 
problems involved was started by the investigation of 
the flow of a hot fluid over a porous flat plate under the 
conditions of uniform fluid injection from the bottom 
The assumptions made in the present 
(1) the inverse proportion between 


of the plate. 
investigation are: 
mass density and temperature inside the boundary layer 
is used, and the viscosity is assumed to be proportional 
to the square root and three-fourths power of the tem- 
perature; (2) the Prandtl Number is assumed to be 
equal to unity; and (3) the flow is assumed to be 
laminar, and the fluid flowing along the wall and the 
coolant flowing through the pores are assumed homog- 


enous. It must be borne in mind that in the turbulent 
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boundary layer that occurs in actual practice, the re- 
duction of heat flow is probably smaller. Moreover, 
the isolating effect of the coolant film is better in the 
case of a laminar boundary layer. Hence, the present 
calculations give an upper limit for the decrease of the 
heat transfer coefficient by sweat-cooling. 

In the solution of the boundary-layer equation, the 
Prandtl Number was assumed to be equal to unity; then 
the temperature 7 could be expressed by a certain 
parabolic function of the velocity only, which satisfied 
both the equation of motion and the energy equation. 
The velocity profile was assuined as a polynomial of the 
fourth degree. 

Knowing that the quantity of heat per unit area re 
moved by the plate is equal to the quantity of heat per 
unit area absorbed by the coolant at the surface of the 
plate, a relation between the wall temperature and the 
amount of coolant needed is established. For a pre 
designated wall temperature, the amount of coolant re 
quired per unit time can be determined provided the 
temperature of the outer boundary layer and of the 
coolant are known. A comparison of the above theo 
retical results with some available experimental data 


shows only fair qualitative agreement. 


THE MOMENTUM AND ENERGY EQUATIONS FOR THI 
BOUNDARY LAYER 


The equation of motion in the boundary layer of a 


flat plate in steady compressible flow is 


Ou Ou O Ou 
pu + pv = vi ) 
On Ov Ov\ Oy 


where the x-axis is taken along the plate in the direc 
tion of the free stream, the y-axis is perpendicular to 
the plate, and w and v indicate the x and y components 
of the velocity at any point. 

The equation of continuity is 


(O/Ox) (pu) + (O Oy) (pv) 0) 2) 


Che equation that determines the energy balance be 
tween the heat produced by viscous dissipation and the 
heat transferred by conduction and convection is given 
in the following form: 

0 = O ; O or On \- 
pu (cpl) + pv 8 (« ) t u( ) (3 

Ox Ov Ov\ Ov Oy 
In reference |, it has been shown that, for a Prandtl 
Number equal to unity, Eqs. (1) and (3) can be satis 
fied by equating the temperature 7° to a certain para 
bolic funetion of the velocity “ only. This relation 


between 7’ and 1 is 


F ie Dw u k | u u 

“tho (. \) t wey ) } 
rT i 7 l 2 l U 

The perfect gas equation leads to (for zero pressure 


gradient): 
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pon = 1,7 5 

It will be assumed that the variation of the viscosity 
with temperature is given by 

wile T/T)? 6 


For air at normal temperature, the exponent is about 
0.76 instead of | 2, but it increases slightly with 7). 


Integrating Eq. (1) between the limits y = 0 and 
y = 6, utilizing the equation of continuity, and satis 
fying the boundary conditions, y = 0; u = 0,7 
y = 6; v = v(6), u {’, one obtains: 


3h . na [ ae . on ; 
u-dy IUAY + py,~voU = ‘ 
Ov Jo ait, oxo" si ‘ "Oy. 


where p,, is the mass density of the injected coolant 


SOLUTION OF BOUNDARY LAYER EQUATION 


In order to solve Eq. (7) a polynomial of the fourth de 
gree as an approximation to the velocity profile is as 


sumed: 
an + hn? $ cn” + dn’ 


where 7 = y 6. 
The values of the quantities a, b, c, d can be calcu 
lated from the following boundary conditions: 


U, O7u/dy? 0, Ou/doy = 0 


_ (ou ~~) 
pain (3) 7 »» (5 
Ou (*" 
(5°) a 


Chis gives 


2 OX 2-+ SA | iA 
a D i P d= 
I+A I1+A I+A | \ 
)) 
V6 ‘y” . . . 
and \ : T,. Hence, the velocity distribution 


Oy 
in the boundary layer is given by 


2 Gd P+ 8... 14+ 3d 


The density p and viscosity uw are expressed as 10! 


lows: 
i i k | 
p a (= I )n + 5 Ven | 
ind 
i. ie k— | — a 
; mE — (50 -1)n $5 ea »| 


where ” 1/. or 3/4. 
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5 Upon substitution of Eqs. (10), (11), and (12) into the first two terms of Eq. (7), one obtains: 
si ( : [ | —_— | @ e + 68 + Gy) : T, (2 ’ + 68 + By) X 
_ “uay = sa—- 27 DO ry oe | Ze Ye + 66 + O 
iro oY = Of — 19/2] 2 "Him. oS ¥ 
6 l (30 + 88 . 3 l dd 13) 
+ (Sa SB + Oy) = ¢¢ + 2 7 20D) (1é 
(1 +)? " “a4 . i+hide * 
about 
: where 7 
() ind poe aa . var nT 6 or 
7) T, Tl, + 27, -— VT,? + 47,4 
Satis € log — + log — 
2 Ty VT72+4+ 47, T, + 274+ VT? + 472 
P Ta 1, + 2Ta — VT,? + 47, 
. = L+ Tre + > = log — 3 aa 
" VT,? + 47 Ty, + 272 + VT)? + 47, 
b3 — (1/2) — Tov on Tae, a (1/3) 7,8 — 7 ay 
| . ae) pind Ma, OS Eee 
T,= 1 + (T, ata Tos) 7 r= V/ " T = / (7, = MT? 
and 
h de 
d ‘ OW, l ) i 3 - 
a pu? dy = 4(12D + GSE + 156F — 180G + 104/7 — 241) x, 
dx Jo “ — Of(k — 1)/2] Wek (1+ )3 
14D + 20E + 36F — 20G — 16/7 + 241 — SJ + T,(12D + GSE + 156F — 180G + 104H — 24])| X 
Ss 
) l 
~+ T; (4D + 20K + 36F — 20G — 16H + 247 ~ SJ) _ + (9D + 48E + 100F 
1 + vA)? (1 + A)§ 
le , ’ . : A , 
ileu 96G — 36H) + (6D + 28E + 44F — 12G — 32H + 24) + (D+ 4E+ 4F4 
(1 + \)? (1 + A)? 
: | \ dr 
1G — SH + 4J/) ; (14) 
(1 + v)2h de 
where 
ae 1 1 T, 
I] a, I[=—-p J | log — + ( + Ty 
2 ie Z 
l T, T rT \ T,? 2 ly l a . 
\ log * G + Tha Tap 
iN VT,? + 4Ta T, + 27a + VTi + 470 f 
‘ Be aaa ais Ett . te . —_ 
Q) I _+7%G+ Taya, E = Pp .F —-TG, D=-.—-T,E+ Taf 
~ ’ ‘ 
_ Hence, Eq. (7) becomes. 
Gp l 2 . N l ae *. N 
= A, t (B, t Tl, Ay) ee ae TB, r Cy = ine dD; : he 
O[(k 1) /2]M:? (1 + A)? (1 + A)’ (1 + A)* (1 + A) (1 + A) 
i() . ee l : r : l ‘|dd ; l P 
(hy —_ T ,H;) Fy — Gy Ol l t pons = 0) (15) 
(1 +A (1 + A) 1 + AJdx 3(4 + 1) (A + Ty,) 
Ol Introducing the nondimensional length & (QO/1V)*(Ux/), Eq. (15) is then reduced to 
dé LS 7 (A + T,) ae nN 
} = | = 4 T (B, T TA, ‘ 
ah [(k — 1)/2|M?py, L8(A + 1) (A+ T,) +1 (1+ \X (1 + X) 
| d* rN l 
7B, T Ci T dD, T (fF Tt) Fr G| (16) 
(1 + A)? (1 + A) (1 +A 1+A 
? wher 
A 12D + GSE + 156F 180G + 104/77 — 24] 
B 1D + 20E + 36F 20G 1677 + 247 + SJ 
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C, = 9D + 48E + 100F — 96G + 36H 
D, = 6D + 28E + 44F — 12G — 32H + 24] — 2a + 2 — 68 — 64 
Eki = D+4E+4F + 4G — 8H + 4/ 


F, = 3a + 8B + 67 
Gi; = a+ 26 + 2 
H, = 2a + 68 — 2e + 6y 


Since Eq. (16) is linear, it can be solved by direct integration. The solution is 


LS » {MM P, 1—T; 
¢ a kh i — } £1341 — T)? Pi + Qi — T)R:] X 
(ek — 1)/2)M? wy '3 lL +X 2 
bby | a Wy — ar ; 
log =; — = + V3 (v, + (1 — 7)? (877 — 67, + DP, 
3X2 + 3(1 + T,/)A + (387; + 1) ‘3s 3 
/3}2 + (1 T;)} 
(1 — 7)/)?Q, — 7,/3 (877 — 6T; 4 DR, | tan ¥ >! = + ¢ ¥ 
V1 + 67, — 37/7 V1 + 67, — 3T/? 
where 
Ny | nrk l ar Ny | nrk | 
T, t M2(p~v/T), T;= 7 M2? (up~ TT" 
12 12 2 S 2 
M, C; F, = A, t dD, . C: : G 
N, = B, 27 2A, Tt ky - Ci D;,, P, BB, A, 
—O; = 3A, 2B; Ki C; t D,, R, = A, IT, 
Che constant of integration C’ determined from the initial condition & = U, \ = T, is thus, 
-18 ” : M ; l— T,) A ia 
"ag (* )} Tr, Pn [,) FP, + = 194 T,)?P; + QO: — TyRi| X 
((k — 1)/2).W?\u,/' 3 2 ' 
2 M 
log rs) \ LQ “*) + (I L',)*(31 7 l) P - T;)? Q 
i | 301 — T, | 
; air 7’, + DR: | p tan! : { 
3 V1+ 67, 317 V1+ 67, — 7 


It may be noted that Eq. (17) is not valid for a Mach Number equal to zero. However, a modified equatior 
corresponding to Eq. (17) for Mach Number equal to zero can be derived as follows: 


For a Mach Number equal to zero, the differential equation in £ and \ corresponding to Eq, | 16) is reduced ti 


dé 1S Mir A + iy q nN a nN 
- 7. ( : A + (By + 7T,A2) — + 
dn Tins sa + 1) a+ 7) + 1 i+ (1 +): 
B nN \ | 
1 ( D / Pst } (7 Is 
1+A7 l d) | \ \ 
where 
Ag = 12d. + 68h. + 1562 LSOd. + 1046 D4 fo 
Bs = 4ag + 20d. + 36c2 — 20d, — Ibe. + 24f Ng2 
Cy = 9a. + 48h. + 1006. — 96d, + 36e. 
Dy = 6d. + 28h. + 440 12d. — 32e. + 24f Pao + Ve OD» Lae 
Ey = dg + 4b. + 4a. + 4de Se. + 4go 
Fs = 38ag + 88: + bye 
Go = a + 282 + 2e 


IT, = 2a 2e + 68 t by 


a 
yo 
r 
~ 
' 
~ 
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fe ao ; fe 39 - zs ? 
y es >i ae Ps 
i. | 1 
— a, * ad 
) T; - = 8) 1 - 7 
| Te I l. , 
4 = e ( ( T De 
7 T, — T, S T; 1 


Phe linear differential equation [Eq. 


1S)] ean be solved in a manner similar to that of Eq. 


16), and the solu- 


tion 1s 
ISn, (“\)\" P, (l— 7 a. 
t = ~ d i + (1 — 7,)? Po + Qe — T,R2 
o ae ~ Tats i+ 2 ; ‘ Ral X 
! Lt. + V31 -(.N vs) 1 — T,)? (37, lr, +1)P 
og . va na T t LVe — | - *(37,° — 61, + 2 
32+ 3(1 + T,)A + (87, + 1 Li 3 carne we ie 
ates Ris . l V3[2 + (1 + T,) 
(1 — T,)?Q@ — — (3T,? — 67, + 1) Re| 7 tan-' — ol i + C” (19) 
3 V1+4+ 67, — 37, V1+4+ 6T, — 37, 
where 
@ (my — 1)/1l2tforp~ VT 
T, = (mr — 1)/8 for np ~ T” 
VW. = Fe— Co = C.+ G. — A, — Dy, Ne = 24.+ D7, - B,- hk, - C 


P, A» - Bo, Q2 = 2B. 4 


rhe constant of integration C” in Eq. (19) is determined from the initial condition § = 0, \ = 


18 d 10 Me . r .~2Z 


ty — Mi 


val 5(: ie) + (1 


rhe solutions as calculated from Eqs. (17) and (19) 
are shown in Fig. 1. It was noted that the relation 
between £ and ), is linear except in the region — < 1.0. 
On examining Eqs. (17) and (19), it is evident that the 
first term on the right-hand side of the equations is 
linear in A, and this dominates the solution when & > 
1.0. This linear relationship between \ and & can be 
shown by their constant slope when > 1.0. The 
growth in boundary-layer thickness with the increase 
of Mach Number and temperature ratio (outer bound- 
iry-layer temperature to wall temperature) are shown 
in Figs. 1 and 2. Since the effect of compressibility is 
to increase the heat transfer through the wall and since 
the amount of heat produced in the boundary layer 
increases with speed, the effects of both the increase of 
Mach Number and temperature ratio to the boundary- 
laver thickness are the same. 


In Figs. 3 and 4, velocity profiles were plotted against 
y(v vy) and yU/y, at various — and Mach Numbers, 
respectively. Because of the injected velocity of the 
coolant along the plate, inflection points appeared in 


T,)(8T,? — 6T, + 1)P2 + (1 


C, + E, — 3A, — Do, Re = He — Az 


—T,, as 


* (3(1 — T,)? Ps + Qo — T,R2] log (1 — T,)? + 


7 
- T,)°Q2 — — (3T,? — 6T, + DR| x 


0 
V3(1 — T,) _ 
V1 + 67, — 87, 


| 
tan~! 


V1 + 67, — 3T,? 
the velocity profiles. It was found that the larger the 
value of v/U, the farther the inflection point moves 
outward from the plate. The instability of the laminar 
boundary layer may be interpreted from the inflection 
points occurring in these curves. 

The temperature distribution at different v/U ratios 
and Mach numbers are shown in Figs. 5 and 6. It is 
noted that the temperature gradient at the wall in- 
creases as the Mach Number increases and decreases as 
vo/U increases. This behavior indicates that the heat 
transfer through the wall increases as the compressi 
bility of the flow becomes more pronounced and de- 
creases as the injection of coolant increases. 

In Fig. 7, a comparison of velocity profiles between 
the present method and the Karm4n-Tsien method,’ 
based on the similarity law, was made. The results 
show good agreement, although the present method 
gives slightly less heat transfer at the wall. 

The above calculations are based on the assumption 
that the viscosity u is proportional to the square root 
of the temperature; however, no appreciable differences 
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BOUNDARY LAYER ALONG A POROUS FLAT PLATE 
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in the results are noted when the viscosity is assumed — or 
to vary with the three-fourths power of the temperature 
ratio, except for large Mach Number. m _ tin ~ i (% ) 
T; — ls. Cr D4 
\\ \LL TEMPERATURE AND COOLANT DISCHARGE tem I ions 
nr(l + : 1/7?) l 


fhe temperature field in the laminar boundary layer - F 
ax . 
o (A+ 


described in the previous section may be used to deter- 
mine the amount of coolant required to cool the wall to 1) (A + Ty) 
a predesignated temperature from the relation stating 
that the amount of heat per unit area removed by the AS - 
plate is equal to the heat absorbed by the coolant. MN — Te - a(t) | f dx 

T To ¢ Ji Jo (A+ 1)(A + Ty) 


Thus, one obtains 
EET aide. af , 
» * tae dx = pyovoty (L — To)k (<0) for lJ = 0, where 7) is the temperature of the coolant 
” e u -y ‘y . . - - 
and nr = 7;/T,,. On the substitution of \ for x from 
Eqs. (16) and (18) in the above integral, one obtains 


for \J > O, and 


(22) 


[ dx | Spin (ey 2] 1 [° | | - 
Jo XY+1 - Ww] » J -T7 3d? + 3(1 + T/A + (1 + 3BT)) 


(AX + Ty) [((k — 1)/2] AP yy 
l 

M (VN, + 7;R1) = me EE T;) ja 23) 

jan + 1 v1 +»? 7 a+ypsf ( 


for \J > 0, and 
[ dx Snr be (ey “| f l | 
Jo A’tDOAHT TZ hae = Dw W/ », To_ 3X2 + 3(1 + T,)A + (1 + 37,) 
7 _ l ” P, 
M, + (No + T,Re2) -— (1 — 7,) - -ldd (24) 
(1 + A)? (1 + A)# 


for JJ = 0. 
Substituting the above integral into Eqs. (21) and (22), respectively, the following expressions for the tempera- 


for \J > 0, 


ture ratios are obtained: 


uw \fk — | 
“p MW “ 
iy — Ze (« Ct )( 2 ru my l) by QO 2 — a 
i iy T —=™ k _ ( S (Is0 - T,)*P, - (N, 4 T,R;) | a . +- = - x 
- 0 of mo( T r wr) _ | Kur 9 


P, 3, (1 + A)? 
(1 —- 7,)(Ni + 7,R:) + 1 — Ty) [1 —3(1 — T,)?] Pi} log + 
a+»? si 9! f 1 1) f | 3 (1 7°;) | P\} log 3x2 + 311 — T)dA + (37, + D T 
, 2 
V3) M, + (83T/ — 6T, + 1) (Ni + T,Ri) + 3801 — T;))? [1 — (87 — 6T; + W1Pip x 
” 


| V3f2A4+ (+ zal)" | ! . 
. tan! (25) 
V1+ 67, -— 3Tf V1+ 67; — 3Tf Ts 


and for M = 0, 


Mr 
eo — Ri(nr — 1 
a Ty —o (")(Sy [301 — T,)?P, — (N T,R 3 b Sink! } 
i To 6n7 Mod \W, e pts mae ee a + 2 (1 + k)* 
(1 + A)? 
t 


J 
2\ 


(1 — 7,) (Ne + 7,R2) + (1 — T,) [1 —3(1 — T,)? Pin y 
, ie oI (1 — Ty)"] Paplog 2 43 (1 + T,)s + (37, + 1) 
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The limit A, can be determined from the curves, A vs. & 


for a corresponding £,, provided the Reynolds Number 
and Q/W are given. 

In Figs. 8, 9, and 10, the temperature ratio (7; 
To) /(T 
flow ratio at different Mach Numbers, Reynolds Num 


Ty) was plotted against the coolant mass 
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bers, and viscosity variation with temperature. For 
predesignated wall temperature, it is shown that the 
rate of coolant discharge increases as the Mach Num 
ber increases and decreases as the Reynolds Number 
increases. It is also shown that larger amount of cool 
ant is required to cool the wall to a predesignated tem 
o rytl/eo - ryt / r 
perature foru ~ 7’*thanforu ~ T’*. Theexperimen 
tal data, taken from reference 3, were compared with 
the present calculations (see Fig. 9). The discrepancy 
between the theoretical calculation and the experi 
mental results is mainly due to the fact that the experi 
mental conditions do not correspond to those assumnied 
in the theories. 
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Shear Distribution in Beams with 
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SUMMARY 


The general equations for the shear flows in horizontal panels 
of a horizontal beam with m + 1 stringers are indicated. The 
three-panel beam with constant stringer areas is solved in detail, 
und the results are shown by curves in terms of the parameters p 
ratio of outer web thicknesses to center web thickness) and b (> = 
vt Ah with 1 = beam length, ¢ = outer web thickness, A = 
stringer area, and kh = panel height). Tests of a three-panel 
cantilever beam with a concentrated load verify the calculated 
results insofar as the test beam represented the assumptions used 


in the theory. 


INTRODUCTION 


ONSIDER A CANTILEVER BEAM (Fig. |) that carries a 
distributed load w(x) Ibs. per unit length. Let the 
horizontal webs ¢;(x) (4 = 1, 2,...., 2) which connect 
the horizontal stringers A,(x) of the beam be thin, so 
that the webs take only shear stresses with a uniform 
distribution in each web on a vertical cross section and 
the stringers take only normal stresses. Let the webs 
act in diagonal tension, and assume a sufficient number 
of vertical stiffeners to be present to support the buckled 
webs and to make any necessary redistribution of shear 
loads between the horizontal panels. Assume the 
heights h; of the horizontal panels to be constant, and 
let gi(x) be the shear flow in pounds per inch in the 
panels. 

Ordinarily, in the calculations of the stringer normal 
stresses and the web shear stresses for the beam of Fig. 
1, the bending and shear deflections of the beam are 
neglected and the usual normal stress formula s = 
My/I is used. However, in many aircraft beams, es- 
pecially those with cutouts, the equivalent thickness of 
the webs may vary widely from panel to panel so that 
deflections cannot be neglected in the distribution of the 
loads. Such a beam is the side of a fuselage that has a 
row of windows with light supporting frames; the 
frames can be replaced by a thin web of equivalent stiff- 
ness, and the analysis carried out on the resulting con- 
tinuous structure. Local effects of the cutouts can be 
analyzed separately and combined with the general 
solution. With suitable modifications, the results are 
applicable to beams with supports other than cantilever. 


DERIVATION OF EQUATIONS FOR THE SHEAR FLow 
From the well-known relations between stresses and 
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displacements, the following relations can be written 
(¢ mw 1,2 ..%.., 8) 


. 


S$; = E; (du;,/dx) (1/A;) / (q;-1 — Qq@dx (1) 


Ga | owt) 1 “ my 
w= 641 (5) + (5. 


where s; = the normal stress in stringer 7, u; = hori- 

zontal displacement of stringer 1, w; = horizontal dis- 

placement of web 7, v, = vertical displacement of web 

i, Eis Young's modulus, and G is the shear modulus. 
From the continuity of the structure, 


Ow; 
Oy 


ui)/hi, Vi = Y 
and Eq. (2) can be written 
Upo41— U2, — Uy; 


Si - Sig: — 
q h i a hy 


withz = 2,3,...., m, where S; = 1/Gjt; Differentiate 
Eq. (3) with respect to x and substitute from Eq. (1) to 
get 
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l 
(Sigi) — B; f gdx + pv, f gi —; dx + 
ax a 


d 
F; faa ,,dx = (Sigs) — B, / qdx + - (4) 
F dx 


F, / Gia, £= 2,3, 22% x nN 


where 


h.B, = (1/AE;) + (1/4; 414,41), Di = 1/hAiE,, 


F, = 1/h Ajo Es +1 


The n shear flows g(a = 1, 2,...., m) can be deter 
mined from the ” - 


equilibrium condition 


> % hq — 
s=1 


/ wx)dx = P (5 


« 


According to ordinary differential equation theory, the 


set of m equations in Eqs. (4) and (5) may be expanded 
into 2(m — 1) first-order linear equations plus two zero 
order equations, or the set may be reduced to a 2(n — 1) 
order linear equation for one of the variables. For the 
general case set up here, the equations have variable 
coefficients so that they are, in general, difficult to 
solve. It is considered more satisfactory to start directly 
from Eqs. (4) and (5) for each particular type of prob- 
lem. 

The constants of integration may be determined from 
the boundary conditions on the beam for the stresses 


and displacements at the ends. 


BEAMS WITH CONSTANT DIMENSIONS 
For the general case of constant A,, f;, &;, Gi, Eq. (4) 
inay be differentiated a second time with respect to x 


and written 


(S:D? — Bia + Figz (S,D? Biq — Dig; — 1 - 
Fiqi41 = 0 (6) 
t= 2,3,...., %, where D? = d*/dx*. The set of alge 


braic equations [Eqs. (5) and (6)] may be solved 
directly for the g; by determinants or otherwise with the 
usual interpretation for the operator D. 

For the special case of two panels, Eqs. (5) and (6) 


give 
qe (P — hyq)/he, (D* K*)q, = (K2?D? —K,*) P 
where 

MAK? = By + Dot h(i + Bo) /hy 

MK.? = So/ho, MK,? = (F,; + Be)/he 

M = ay + (hy S» hy) 


This result for the two-panel case agrees with that on 
page 9 of reference 1, provided vertical webs are used 


there. 


THE THREE-PANEL BEAM 


Consider the three-panel case for which Eqs. (5) and 


AERONAUTICAL 


1 equations in Eq. (4) and the 
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6) give 
hsqz; = P — Inq, — hege 
(Ko — Ada = — [(Ai + Ks)D? - | ‘ 


K. + Keg) \q2 + [K7D? 


- (K3 + Ks)] P| 
D' L,D? + Le) qe = (L,D4 - 


L4D? + L;)P 
where 


3S, Ko = h3(B, + Dz) — hyF» 

h3S;\ Ke. = hg(F, + Bo) + heFp 
hyS,K; = Fo, haJKg = hgB, + h,Bz, 
hsJ = h3S; + S3, hgJKs = heS; 
haJ Ke = h3(D3 — Fi) + A2B; 

hJK; = Ss, hsJKs = B; 

K, + A;)Li = K,iK4 + AoKs + Ko + Ke 

(K, + K;)L2 = KoKg + KoKes, (KA, + A;5)L3 = K: 
K, + K5)L4 = AoKi + Kz +Kg 

AK, + Ks)L5 = KoKs + Az Kg 


SiKi -- So 


From Eq. (7), for P, a concentrated load at a distance 
a from the fixed end and a distance / from the free end 
l= a+ d) of the beam, 





Jou = Cig sinh mx + Co, sinh mex + 
C3, cosh myx + Cy, cosh mex + (L;P/Le) . 
fia(mix) + fog(mex) + (L5P — 
2b = tiv (9114x ) + fos (Mox ) 
where 
Im? = L; — (-—1)' VL — 4b, 2 ee 
Put Eq. (8) into Eq. (7) to get 
dia = JSilmy)f a(x) + Ji(mes)fog(mex) + JoP ) 9 
° » ie 
di = JSi(m)f w(x) + Ji (me) fo,(mox) 
dza = HWy(m)fiq(mix) + Hy(mz)fog(mox) + HeP | (10) | 
dao = E(m)f (mx) + Hime) fo, (mex) 
where 
(Ky — Ky)Ji(m;) = [(Ay 4+ Km? ~— (K, + Ks), 
t=i12 
Li(Ko — K4a)Jo = Ls(Ke + Ao) Li(Ks + Ks) 


hsfl,(m;) = —he — hy, (m;) 
h;H> =]|—- h, J» - (hols Lo) 


The case of A;, /,, ti, hi, G,, having different constant 
values in region a and 6 can be solved to give similar re- 
sults. 

CONSTANTS OF INTEGRATION 


The boundary conditions u;, = Oat x = 0, s;, = Oat 


x = 1, u,, = uy at = a, and s;,, = Sy atx = 6 
(1 = 1, 2, 3, 4) give, from Eqs. (1) and (3), 
Sidia = S2Gaqa = S3G3q, atx = 0 
Siap' = Sedu’ = Say’, atx =! 
Sildia — Qin) = S2(Goa — Gov) = S3(G3a — Q30)> 11) 
atx=a 
Sildia’ — Qi’) = S2(doa’ — Goo’) = S3(qza’ — 
@y,); atx=ea 
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Put Eqs. (8), (9), and (10) into Eq. (11) to determine 


MCz, Mi(m2)P, MC,, = M,(m,)P 
Cap = C;, (1 — cosh ma) 
Ci = — C,,(1 cosh ma) tanh m,/ | 
ce. C3, [1 — cosh ma) tanh m/ + 
sinh m2,< | (12) 
Cu = Cy, (1 — cosh ma) 
Cop - — C,,(1 — cosh mea) tanh m,/ 


Co = — Cya [CU — cosh moa) tanh ml + 
sinh mya] | 
where 


Si,Ji(me2) — Se 
S3fT\ (mz) — Se 


Si) J1(m) - S» 
S3f7,(m,) — S» 


SiJi(m;) — Sz 
S;H\(m;) — S» 


Sols — S,L2J. 


on = ee She 


SPECIAL VALUES FOR BEAM CONSTANTS 


As a particular case, let 


h=h, = hs = 2h,/3 
Asda hen hw hk ; 
Satan min hake head 
t= th = ts = phG/G = plrg 


where p allows for f, to be different thickness from 4, and 
t;and also may allow for a reduced G2 caused by diagonal 
tension in fo. 

For this case Eqs. (8), (9), and (10) can be written in 


the form 


107hqj. : 5 
p = N;,R,(bkix /1) NioR,(bkox /l) + Q; 

LOThgi» A , (14) 
p = N,,Rp(bkyx/l) — NjoR,(bkex/1) 

tau 3.2.3 

where 


Q, = 26, Qe = 36, Q3 = 30, b = U(t/Ah) 

(5p + 2)k2 = (13p + 37) — (—1)'T 

I? = (13p + 37)? — 321(5p + 2) 

5p + 2)T No; = 1170p? — 3371p + 1,867 + (90p — 

71)(—1)'T 

Ni; = [(20/3) —-(5p + 2)k? | No; 

N3; = —[(46/9) — (2/3)(5p + 2)k?1]No; 

R,(bkix/l) = cosh (bkix/l) — U; sinh (bkix/1) 

R,(bkix/l) = [1 — cosh (bk,a/l)|[cosh (bkix/l) — 
tanh bk; sinh (bk,x/1) | 

U; = [1 — cosh (bk,a/1)| tanh bk; + sinh (bk,a/1) 


THE PARAMETERS Pp AND 0 


Eq. (14) presents the shear flows in terms of the two 
parameters p and b = I(t/Ah)'*. It is noted that a 
similar result holds for any initial proportions between 
the constants /;, A;, E;, G; andt; Since p is a measure 
of the diagonal tension and of the relative thickness of 
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the f2 web and since 6 has the length / in it, as well as 
t, A, and h, the parameters p and 6 are considered to 
represent the most important effects that the deforma- 
tion of the structure has upon the distribution of the 
loads. Note that the Q; in Eq. (14) are the shear flows 
when all deformation is neglected. 

Fig. 2 shows the effects of the parameter p upon the 
shear distribution for a selected value of b and a/I = 1, 
while Fig. 3 shows the effect of the parameter 6 for a 
selected value of panda// = 1. Fig. 4 shows the shear 
distribution for the load located at different points on 
the beam. Note that the shear flows are not zero be 
tween the load and the free end. 
comparison of the shear flows in this case to the case 


Fig. 5 shows the 


when the shear flows are taken as zero outboard of the 


load. 


Tue Test BEAM* 


In order to compare the theory for a three-panel beam 
to an actual beam, a test beam, with the beam constants 
related as in the above sections, was built. For this 
test beam, h = 10 in., A = 0.62 sq.in., EF = 3G = 10° 
Ibs. per sq.in., £ = 0.081 in., fe = 0.020 in., / = 100 in., 
the material being 24ST extruded angles for the string 
ers, 24ST Alclad webs, and 24ST extruded angle stiff 
Without diagonal tension the 
Fig. 6 shows the 


eners spaced at 8.33 in. 
factor p = 4, and the factor 6 = 11.43. 
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Effect of parameter "pon shear 
Fig. 2 


* The author wishes to acknowledge the aid of Asst. Prof. V. A 
Valey, of the United States Air Force Institute of Technology, in 


designing and testing the test beam. 
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test values obtained by the use of SR-4 rosette strain 
gages as compared to the calculated values. It is noted 
that g; is low at the free end and high at the wall, while 
g; is high at free end and low at the wali. This effect 
was due to the method of loading and the supporting 
jig. The load was applied at the free end as a concen 
trated tension load into a light end stiffener (Fig. 2), 
which tended to distribute the load to panel 38 first. 
The load was taken out at the wall through a light 
stiffener in tension, which tended to concentrate the 
load in panel 1. ge is high toward the free end because 
of the fact that the assumption that the webs take only 
shear is not obeyed in this test specimen. The webs 
if, = ts = 0.081 have areas equivalent in magnitude to 
the stringer area, so that at low bending stresses near the 
free end the webs still act as stringer areas, which will 
necessarily raise the shear stress in the middle web /:. 
rhe strain gages indicated that the /; web carried com 
pression stress according to J/y// in the outer half of 
the beam but that this stress was below J/y/J in the 
inner half of the beam and actually decreased from its 
peak value toward the wall. The ¢; web carried tension 
stresses approximately according to .Wy/J throughout 


its length. 


CONSIDERATIONS ON APPLICABILITY TO ACTUAI 
AIRCRAFT STRUCTURES 


In design and stress analysis of an aircraft structure 
with variable webs, the three-panel beam analysis is 
considered to be sufficiently good for most cases. Re- 
gardless of the number of stringers, the structure can be 
regarded as made up of three representative sections for 


the middle, the top, and the bottom. Once the shear 
flows are obtained for the representative three-panel 
structure, the analysis may proceed by superposition of 
the usual section property results — the shear flow in the 
middle web represents a certain per cent of the bending 
moment as carried by the beam as an entire unit, while 
the gq; and qs; shears represent the per cent moment 
carried by the beams with web | and 3. The final loads 
in the various members are obtained by superposition 
of the results for the three representative beams. De 
termine the parameters / and 4 for the structure, and, 
from curves similar to Figs. 2 and 3 for the type of load 
ing on the structure, divide the bending moment into 
components for each beam. Calculate the loads in each 


beam, and add the results. 


CONCLUSIONS 


The use of the usual \/y,/ formulas for aircraft 
beams with variable webs may lead to misleading results 
in many cases. The deformation of the structure and 
the consequent redistribution of loads can be repre 
sented in terms of two parameters, P and 6 = IV t/A h, by 
families of curves. By judicious use of such curves, a 
more correct analysis of the structure can be made, and 
it is believed that lighter and better proportioned 


structures can be designed. 
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spondents 
—— 
Note on Supersonic Wing Pressure Distribution 


Ellis Lapin 
Douglas Aircraft Company, Inc., Santa Monica, Calif. 
September 6, 1949 


[' THE PRESSURE DIFFERENCE at any point p(é,,y,) of a flat 
wing of trapezoidal half-plan form with sweptforward trailing 


edges is denoted by 


B 2 Bu Bua E BAun B Aou 
* 2. SS —— a r = + : (1 
la T a l a I a I a 
where 
Bus ie me 
= pressure difference produced by an infinite delta 
Ca 
wing 
BAyu ; ; . 
7 = incremental pressure difference induced by the pres- 
a 
ence of side edges 
8 Acu ‘ 
7 = incremental pressure difference induced by the 
a 
presence of subsonic rear edges 
8? = M?—-1 
nt = coordinates oriented along Mach waves from wing 
vertex 
a = angle of attack, rad. 


then by application of the method described by Goodman! it may 
be shown that the incremental pressure difference, 8B Asu/ Ua, in 
the regions adjacent to the subsonic edges, as illustrated in Fig. 1, 
is given by: 


“hy 





~“ 


rm. (2 A 1 
7's, P = _ Bua dn 
Ue Ua (2a) 
"se 


(4#') _ Vw— ali 
Ua 3 T AE —— 
- he nV ar(é,) = 


(BAou/Ua); = (8 Anu/Ua)s (2b) 


(22 (2) ea) \ ‘~ — nr(to) . 
Ualyw \ Ua Ua /s 7 
“np (é : 
/ (8 Aou/Ua)s dn (Qe 
/7 
” % — 2) V nr(é) =-9 
(22) (A=) V 10 — nr(to) 
Va 17 la I r 
“ng (& . 
/ (3 Aou/ Ua) i3 dy 2d) 
Jy] —S V nr(é) a 
(ae (es) ve - aft 
Va 29 in la ? 7 
“ap (E : 
(8 Aou/ UVe)is dy Qe) 
n 
vill % — ”) Vv nr (tp) — 7 


etc., where np(&) is the point on the trailing edge with & = & and 
ns(&) is the point on the extension of the side edge with — = §). 
It can be demonstrated that the pressure difference, as computed 
by means of Eqs. (1) and (2), complies with the Kutta condition, 
AC, = Oat the trailing edge. In the regions more remote from, 
but still influenced by, the trailing edges, the incremental pressure 
difference, 8 Aou/ Ua, may be found by application of a superposi 

tion principle that accounts for the pressure resulting from the 


Mach wave reflections: 


(BAwu/Ua)e = (BAU, Ua 3a) 
(BAwu/Ua)yo = (BAu/Ua)s + (BAu/Ua)s (3b) 
(BAou/Ua)y, = (8 Au Ua); + (BAsu/Ua)s (3c) 
(BAou/Ua), = (BAu/ Ualy + (BA u/ Ua): (3d) 


(? Aou B st) (° =) . 
tt Ge Fe Ua js Va I ces 
(fae * (Pas) (2) | 
Ua Jo Ua Jrs Ua Js 
(° Aou ) B Aou | 
= — oa 3f) 
| Ua 18 ( UVa 3 
etc. 


The symbol (8 A,u/Ua), is to be construed as meaning “‘com 


( pb ae) 
Ua 16 


pute the incremental pressure difference (8 A.u/Ua) at any point 
p(n) as though the point lie in region v.”’ 

Since the “‘basic’’ pressure difference, 8ua/Ua, and the side 
edge induced incremental pressure difference, 8 A,u/Ua, repre- 
sent known quantities (cf., reference 2) the computation of the 
pressure at any point influenced by the trailing edge involves the 
evaluation of Eq. (2a), the evaluation of interaction terms such 
as the third term in Eq. (2c), and the application of a special super 
position principle. The form of the equations appears suitable 
for evaluation by means of automatic computing machinery. 
The restrictions of linearized inviscid supersonic flow theory 


apply. 
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Heat Source in a Uniform Flow* 


H. S. Tsient and Milton Beilockt 
Massachusetts Institute of Technology, Cambridge 
August 25, 1949 


Vass. 


_— PROBLEM OF THE disturbances produced and the forces 
acting on a heat source in a compressible fluid has been the 
subject of a recent investigation. The solution obtained is for a 
line source of heat of strength QY (in mechanical units) per unit 
length per unit time in a two-dimensional unbounded flow in the 
x, v plane with the undisturbed velocity Ul’ in the x direction 
The line source is perpendicular to the x, y plane and is located at 
the origin. By neglecting the effects of viscosity and heat con- 
duction and by using the Fourier integral technique after linear- 
izing the partial differential equation, the following results are 


obtained: 


WU = free-stream Mach Number 


p® = free-stream static pressure 
y = ratio of specific heats 
u’ = x-component of disturbance velocity 


= y-component of disturbance velocity 


pressure disturbance, change of static pressure from the 
value p° 
p’ = density disturbance 
a° = velocity of sound in undisturbed flow 


For subsonic flow, M < 1, 
iy = DC x 
© Seep A) — MP xt t+ (1 


,_ v¥-)WV1- are y 


2ryp” x? + (I WP )y 
(y -— 1)MQ x 
? 
f 27a° V1 — M2 x? + (1 — M*)y’ 
(y -— 1)MQ \ + - 10 
a ~ B(w) I(x 
2a’? 4/1 — M2 x? + (1 — My? a? 
where 6(y) is the 6-function defined as 
iy) = 0, vy >0O0 
6(yidy > nee 1 
ind J(x) is the unit-step function defined as 
I(x) = O forx < 0 
I(x) = lforx > 0 
For supersonic flow, > 1, vy > 0 
Y= 1)Q 
yp? / Mt — 1 . 
y¥~ 1@ 
ap be - VM ly 
; (y -— I) MQ 
p 7 = Bx — VV MN? ly 
20° \/ M2 — 1 , 
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ind 


p'(x, — y) = p(x, y 


p'(x, ~ Se — p’(x, ¥y 


Itis then seen that for subsonic free stream the disturbances are 
sourcelike, while for supersonic stream the disturbances are essen 
tially limited to a narrow region along the Mach lines. The 
second terms in the formulas for p’ are, of course, the heat wake 
By taking proper consideration of the linearization, the momen 
tum change due to the heat source is also calculated. This caleu 
lation shows that there is no force acting on the heat source. 

The above results can be immediately applied to the problem 
of arbitrary distribution of heat source by the principle of super 
position. As an example, let there be a distribution of lin 
sources in the x, y plane along a line inclined at the angle @ with 
the x axis. In other words, let the heat source be distributed on a 
plane perpendicular to the x, y plane and making an angle @ with 


the x,z plane. Let the intensity of the source be q per unit area 


perunittime. Then for 1/ < 1, 
iy — Ig sin 
‘ -  -2yp? 1 — MP sin?@ 
y¥— lq cos @ 
yp? 1 — Af? sin? 6 
for Vf > 1, 
; (y lig tan 8 X sin B 
oo Oyp? sin (B — 8) 
; y — lq sin 6 
sii 2yp sin(8 — @ 
where 3 is the Mach angle, sin~! (1/.1/). It is interesting to note 


that for the subsonic free stream, the disturbance velocity vector 
is perpendicular to the source surface, while for supersonic fre: 
stream it is perpendicular to the Mach wave. In both cases, 
the flow is decelerated before reaching the heat source. These 
results have application to the problem of evaluating the true 
flame”"propagation speed from flame front photograph 
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Action of Viscosity on the Supersonic Wave Drag 


E. V. Laitone* 
Special Consultant to Hughes Aircraft Company, California 
september 1, 1949 


SUMMARY 


[' IS SHOWN THAT in the case of two-dimensional supersonic 
flow, viscosity effects can actually decrease the wave drag 
due to thickness because of the finite boundary layer thickness at 


the rear of an airfoil. 


DISCUSSION 


In subsonic flow viscosity effects not only produce a skin fric 
tion drag but also directly produce the so-called pressure drag 
due to the reduction of the pressure recovery at the rear of 
the body. As shown in Fig. 1, there is a flow compression at 
the rear of an airfoil which results in zero drag in nonviscous flow 
but yields a pressure drag whenever viscosity effects produce 4 
boundary layer preventing the stream lines from closing in at the 


trailing edge. 


* Also, Associate Professor, Mechanical Engineering, University of Cali- 


fornia at Berkeley 
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However, in the two-dimensional supersonic flow over an air- 
foil, there is no flow compression, the flow continuing to expand 
over the entire rear portion of the airfoil. The surface pressures 
as shown in Fig. 2 produce the so-called supersonic wave drag due 
to thickness (or pressure drag) even in a nonviscous fluid. Any 
finite boundary-layer thickness would prevent the stream lines 
from closing in at the trailing edge and would therefore reduce the 


nces are pressure drag at the rear end. Therefore, in the case of two-di 
e essen- mensional supersonic flow, the viscosity effects create a drag due 
s. The to skin friction as in the subsonic case but now actually reduce the 
it wake existing supersonic pressure drag, especially if flow separation 


momen takes place. 





Is caleu The two-dimensional supersonic first order lift depends only 
ec. on the mean inclination of the stream lines along the upper and 
oroblem lower surfaces. Since the lift of a two-dimensional supersonic 


f super airfoil does not depend on the flow closing at the trailing edge and 
of line ' the resulting Kutta-Joukowski condition, it is seen that to the 
 O with FIG. | first order the supersonic lift is not affected by viscosity if the 
ed ona boundary-layer thickness is the same on the upper and lower sur- 
> # with a faces. 

Because of this, an improved lift-drag ratio can be obtained for 
an airfoil shaped as in Fig. 3 and having a jet issuing with at least 
free-stream velocity from the trailing edge. If, as in Fig. 3, the 
jet is the same thickness as the airfoil, then the pressure drag 
would be reduced one-half, allowing a higher lift-drag ratio to be 
obtained at lower lift. Generally speaking, this is a special case 
of the Busemann biplane, the total first order lift increment of the 
upper and lower surfaces being the same as that for a single thin 
flat plate of the same plan form. 

The above effects are not encountered for a body of revolution 
in supersonic flow because, as shown in Fig. 4, there is again a 
strong flow compression at the extreme rear of the body of revolu- 
FIG. 2 tion. Then any boundary-layer thickness would increase the 
drag by decreasing this pressure recovery. As seen in Fig. 4, 
the pressure gradient at the rear end is greater in supersonic flow 
than in subsonic flow so that the viscosity effects would produce 
a greater increase in the surface pressure drag in the supersonic 
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An Experimental Verification of Local lsotropy 


Drag 





Stanley Corrsin 

Department of Aeronautics, The Johns Hopkins University, 
Baltimore, Md. 

September 14, 1949 





— MOST SIGNIFICANT IDEA contributed to the problem of 
turbulent shear flow in many years is the hypothesis of 








FIG. 3 
rsonic “local isotropy” due to Kolmogoroff.!) 2, He has suggested that 
» drag — the fine structure in turbulent shear flow may be isotropic. 
ness at Recent experiments of Townsend* * in a wake seem to verify 
‘ this hypothesis. He found that the skewness and flattening 
‘a M<| || factors of the probability density of (aw/at) in the shear flow are 
UX a ! closely equal to those in the (effectively isotropic) turbulence 
wy | far behind a grid. (u = velocity fluctuation in mean flow direc 
n fric ; ‘| tion; ¢ = time.) 
e drag An independent verification of local isotropy has been obtained 
ear of : by the measurement of the shear correlation coefficient as a 
ion at \ \ | function of frequency, 
is flow — 
juce a ‘ . ; Un Un 


~ a | 
at the ’ aie ai 
at the Ain’ Yu’ 


Here, v is the lateral velocity fluctuation; «, and v, are the con 
aici tributions to u(t) and v(t), respectively, in a narrow frequency 


vacsalbd band of nominal frequency » 
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The monotonic decrease to zero of ,Ry, with increasing fre 
quency is a clear indication of the presence of local isotropy 
Fig. 1). The overall shear correlation coefficient, also given 
in the figure, was directly measured by the conventional method. 
Reasonably good agreement with this value is obtained by 


computation from ,Ry, and the power spectrum of uw. 


These measurements were made in the maximum shear region 
of a fully developed round turbulent jet. They are part of a 
broader investigation on turbulent jets® sponsored by the Na- 


tional Advisory Committee for Aeronautics 


Note on Air-Speed Indicator Compressibility 
Error 


Roger D. Schaufele 
Doug las Aircraft Company, Inc., Santa Monica, Calif. 
September 14, 1949 


pw BASIC THEORY of the compressibility error inherent in 
ordinary air-speed indicators has been well established by 
previous investigators.':*}% The present note extends the 
theory into the supersonic region. 

An air-speed indicator is essentially a differential pressure gage 
measuring pitot-static pressure, with a dial marked off in arbi 
trarily chosen speed units. Analysis shows that there is no 
unique relation between the pitot-static pressure difference and 
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true indicated air speed.* However, at a given pressure altitud 
the pitot-static pressure difference is a function only of tru 
indicated air speed, as shown in Fig. 1. An ordinary air-speed 
indicator, therefore, will read true indicated air speed at only on 
pressure altitude. Usually air-speed indicator dials are marked 
off to read true indicated air speed at standard sea-level stati 
pressure. The error in the instrument reading incurred at 
altitude because of the arbitrary use of this sea-level standard 
has become significant in recent years as the performance limita 
tions of aircraft have been extended to higher speeds and al 
titudes. With the advent of supersonic flight, this error becomes 
increasingly important. 

* True Indicated Air Speed (T.I.A.S.) is defined as the true speed times 
also called Equivalent 





the square root of the density ratio | i, Viot/2 


Air Speed (E.A.5S.). 
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Phe analysis is performed for an ideal or ‘‘perfect”’ installation. 
It is assumed that the meter dial has been marked off to read 
true indicated air speed ut sea level, that there is no lag in pres 
sure readings, and that there is no instrument or position error. 
[he pressures are calculated using the compressible fluid rela- 
at supersonic speed, the effect of the bow wave on the 


tionships; 
In addition, it is assumed that the 


total head is included. 
static tap is located in a region of free-stream static pressure 
behind the shock. 

The method of calculation is as follows. 
ititude, the pitot-static pressure difference, corresponding to 


At a given pressure 


a given instrument reading, is converted to a pressure ratio by 
the identity 


Protai/Patatic = [Ptota Potatis Potato) + ] 


From the pressure ratio, it is possible to obtain the Mach Number 
from the relation 
f YO = 2) 
Potal /Petatic = i] 1 + (iy — 3 /2)M?¢ 


for subsonic flight and 


| 7 4:4 (y¥ + D/(y 1 rs 1/(y 1 2¥/(y — 1) 
Wet eee cae ae 


for supersonic flight. Fortunately, tables‘ of these pressure ratios 
is a function of Mach Number have been prepared. Simple 
interpolation will usually give suflicient accuracy. Once the 
Mach Number is determined, the true indicated air speed may 


be found from the equation 


Vig = AV yeRTo 


where 
flight Mach Number 
= ratio of specific heats 
gravitational acceleration 
R = gas constant 
1 absolute temperature 


iltitude density ratio 


The difference between the instrument dial reading and the true 
indicated air speed is the compressibility error, and its negative 
is the compressibility correction increment 

The results of the study are presented in Fig. 2. The com- 
pressibility correction increment is plotted as a function of the 
iir-speed indicator reading. This is the curve that would be 
used by pilots to determine the true indicated air speed from the 


cockpit air-speed indicator reading. Although for present-day 
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AIRSPEED INDICATOR COMPRESSIBILITY CORRECTION CHART 
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aircraft performance the magnitude of the compressibility error 
is not large, perhaps never more than 10 m.p.h., the study shows 
that there is a considerable increase in the error at higher speeds 
and altitudes. For example, an indicator reading of 500 m.p.h. 
at 50,000 ft. would be in error +78 indicated m.p.h. just from 
compressibility effects. 

There are two additional points that arise in the discussion of 
the correction curves. First of all, the corrections are valid 
for any off-standard day so long as pressure altitudes are used 
and the o in the equations is defined as the ratio of the density at 
altitude to sea-level standard density. Secondly, the correction 
curves have a point of inflection at 1J = 1.0. This change is 
caused by the formation of the detached shock wave in front of 
the pitot-static tube and the resulting decrease in total head at 
the pitot tap. At higher Mach Numbers, and consequently 
stronger detached shocks, the compressibility error actually 


decreases. 
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Hadamard’s Solution for an Inclined Body 
in Linearized Supersonic Flow 


E. V. Leitone* 

Associate Professor, Mechanical Engineering, University of 
California at Berkeley 

September 1, 1949 


g gee INCORRECT SOLUTION mentioned in reference 1 for the 
velocity potential for linearized supersonic flow about an in- 
clined slender body of revolution had previously been noticed 
and the correct form given as Eq. (29) of reference 2. This direct 
inalogy to the subsonic doublet was obtained in the supersonic 
case by introducing Hadamard’s concept of the “‘finite part of 
the integral.’’® 

For example, the potential function given by Eq. (11) of refer- 
ence 1 should have been written, in the notation of reference 3, 


* Special consultant to Hughes Aircraft Company. 
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as the finite part of the improper integral resulting from the dif- 


ferentiation so that 


° 3 
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as proved directly in Eq. (8) of reference 2. 
It is also shown by Eq. (34) of reference 2 that 


Ss” (x) ‘ 
g(x) = tal’ S(x 1 + ,? log - + O(r?) 2 
2r x 


where S(x) mr? is the local cross-sectional area of the body of 


revolution. Therefore to the first order (r‘ 0) we have from 


Eqs. (1) and (2) 


, 1 3 
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. 


Note that Eq. (4) contains the first derivative of the cross 
sectional area, while Eq. (3) only contains S and therefore may 
be simpler to use in most cases, especially since it is the direct 
inalog to the subsonic doublet Chis also shows that Eqs. (10 
ind (11) of reference Ll are compared incorrectly 

The corresponding Hadamard solutions for the supersonic 
source distributions for the axial flow along a body of revolution 
» 


are given in Eqs. (16) and (17) of reference 
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Supersonic Airfoil Theory * 


ee S. Wasserman 
Wright-Patterson Air Force Base, Dayton, Ohio 
september 96, 1949 


T' APPEARS that ua better physical picture of Ackeret’s Super 
sonic Theory can be obtained by utilizing the particle velocity 


set up by the plane wave as shown in Fig. 1 
° 
NZ at 
” ¥ 
. 4 © 
> P .90-m+4 <a 
sd are 
= ’ 
piar® 
‘ 
> 90-™ ) 
Fy A 


In order to satisfy the condition normal to the airfoil surfia 
the following equation is applicabk 
V sin p I’ sin (90 f+ 8, 
Che value of U' for a plane sound wave is given by the well-knovw 
equation 


\p »¢ Z 


If @ and pg are assumed small and cos Wis set equal to its valu 


yY 1 C?/V, substitution of Eq. (2) into Eq 1) vield 
Ackeret’s formula 

\p pl \ \/ l } 
* Comments on the Readers’ Forum (tem An Approach to Supersoni 
Airfoil Theory,’ by W. E. Strohmeyer and D. R,. Gero, Journal tt 
Neronautical Scienc Vol 6, No. 8, pp. 506-507, August, 1949 





Combustion Problems in Ram-Jet 


Design 


(Continued from page 713 


If W, = the pounds of air per second passing a unit 
area in the duct and W, = the pounds of fuel injected 


per second, then 


W, = 20W, | frdr 
- H 
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Present Status of Research 


on Boundary-Layer Control 


Continued from page 740 
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